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EXPERIMENTS ON THE EFFECTS IN HETEROZYVGOUS CONDI- 
TION OF SECOND CHROMOSOMES FROM NATURAL POPULA- 
TIONS OF DROSOPHILA WILLISTONI 


By A. R. CorDEIRO 
CoLUMBIA UNIVERSITY, NEW YoRK, AND UNIVERSITY OF PoRtO ALEGRE, BRAZIL 
Communicated by Th. Dobzhansky, March 31, 1952 


The Problem. Natural populations of all species of Drosophila so far 
studied in this respect carry a great amount of concealed genetic varia- 
bility, consisting of recessive mutant genes and chromosomal aberrations 
of the inversion type. Many of the mutant genes are lethal or semilethal 
in homozygous condition.' The effects of these lethals on heterozygous 
flies are insufficiently known, and conflicting views about these effects 
have been expressed by several authors. Berg’ showed that the lethals 
found in natural populations of D. melanogaster are deleterious in hetero- 
zygous condition, and Stern and Novitski* found the same to be true for 
lethals which arise in the laboratory. Dubinin‘ considered, however, the 
natural lethals in the same species to be completely recessive, but he quotes 
the results of Masing who regarded many of the natural lethals to be 
heterotic. In D. pseudoobscura Wright, Dobzhansky and Hovanitz® 
found that autosomal lethals are less frequent in natural populations than 


they would be if they were completely recessive and if the populations 


were wholly panmictic The discrepancy may be due either to deleterious 
effects of most lethals in heterozygotes, or to local inbreeding in natural 
habitats, or to a combination of both causes. Muller® emphasized the 
importance of the viability effects of heterozygosis for deleterious recessive 
genes for an understanding of the dynamics of populations of sexual organ- 
isms, including man. He also referred to unpublished experiments on D. 
melanogaster, which indicated that mutants which are lethal when homo 
zygous are deleterious in heterozygotes. 

The experiments reported in the present article have been carried out 
to investigate the effects on the viability of the flies of heterozygosis for 
different chromosomes found in natural populations of Drosophila willistont. 
Chromosomes which are lethal to homozygotes as well as chromosomes 
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which permit the homozygotes to survive (the ‘““Normals’’) have been 
studied. It appears that combinations of these chromosomes give rise to 
a variety of genotypes which differ in their viabilities under the conditions 
of the experiments. 

Material and Methods. —Samples of wild populations of D. willistont were 
collected at Quarai, Emboaba and Caxias, in the state of Rio Grande do 
Sul, and at Catuni, in the state of Bahia, in Brazil. The effects in homo- 
zygous condition of second chromosomes found in these samples were 
studied as described in a previous publication.’ A total of 52 second 
chromosomes which were lethal in homozygous condition, and 60 second 
chromosomes which produced homozygotes close to normal viability were 
selected. All the lethals were known not to be allelic. 

The lethal and the “Normal” chromosomes are maintained in separate 
strains in balanced condition. A chromosome which contains the domi- 
nant mutant Star, the recessives brown and abbreviated, and a pericentric 
inversion which suppresses crossing over in heterozygotes, is used as a 
balancer.’ 

Four series of experiments were arranged, as shown in table |. In the 
first series, flies which carried Star (and the other genes) in one second 
chromosome, and a lethal in the homologous chromosome, were crossed to 
flies homozygous for a known Normal second chromosome. The offspring 
of this cross should consist of Star flies and wild-type flies (heterozygous 
for the lethal) in approximately equal numbers. In the second series of 
experiments, flies which carried one Normal second, and one Star second, 
chromosome were crossed to flies homozygous for a different Normal 
second chromosome. The offspring of this cross should, theoretically, 
consist also of Star and of wild-type flies in equal numbers. In the third 
series, flies heterozygous for the Star chromosome and for different (non- 
allelic) lethal chromosomes were intercrossed. Since the Star chromosome 
is lethal when homozygous, the offspring of this cross should consist of 
about two-thirds Star and one-third of wild-type flies. The fourth series 
of experiments differed from the third in that the chromosomes opposite 
to Star were known not to be lethal when homozygous. The offspring 
should, theoretically, consist of Star and wild-type flies in a ratio ap- 
proaching 2:1. 

A summary of the results is given in table 1. A total of 359,667 flies 
were counted. They came from 735 cultures, representing 198 different 
combinations of second chromosomes. The cultures were kept at 25°C. 
on banana agar medium; the flies that hatched were counted on alternate 
days until the ninth day after the beginning of the eclosion of adults in a 
given culture. 

Incomplete Recesswity of the Lethals. The theoretical ratios, of 1 Star: 1 


wild-type, or 2 Star: 1 wild-type, expected in the two categories of the 
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crosses can be realized only if the Star and the wild-type flies were equal in 


relative viability under the conditions which obtained in our experiments. 
In reality, the carriers of the Star chromosome were known to be less viable 
than their wild-type sibs.’ 

The ratios of Star: wild-type flies are, however, not uniform in the dif 
ferent crosses (table 1). It can be seen that the wild-type class is less fre 
quent in the offspring of the Star lethal to Normal crosses (55.20 per cent 


TABLE 1 


A SUMMARY OF THE EXPERIMENTAL DATA 
COMBINA MEAN 
TIONS OF NUMBER TOTAL PER CENT FREQ PNCY 
CHROM® OF FLIES WILD OF WILD 
TYPE OF CROSS LOCALITY SOMES CULTURES COUNTED IN TOTAL, TYPE 
Star/lethal X Quarai 32 102 54,724 60 42 60 94 
Normal / Normal Catuni 37 135 100, 202 55.33 54.91 


Emboaba 20) 172 76,500 51 29 55.67 


Total SY 231,516 55.20 

Star/Normal X Quarai j ‘ 979 33 
Normal / Normat Catuni ‘ 24 14,876 16 
Emboaba ; 38 14,716 55.04 


Total P 19,571 57.78 

Star/lethal x Quarai : . 8,251 3.87 
Star/lethal Catuni 2 20,712 26 
Emboaba : 5,068 87 

Caxias : 33,471 10 


Total »? s. 67 ,502 17 
Star/ Normal X Quarai : ) 1,610 
Star / Normal Catuni ¢ 1,760 5.20 
Emboaba j j 3,642 
Caxias : 1 066 


Total 31 ' 11,078 5.32 


in the total) than in the offspring of the Star, Normal to Normal crosses 
(57.78 per cent in the total). The obvious difference between these two 
categories of crosses is that in the first of them the wild-type flies are 
heterozygous for a chromosome which is known to be lethal to homozygotes. 
This suggests that, on the average, the lethal chromosomes used in our 
experiments are deleterious also in heterozygous condition. This sug 
gestion is confirmed by comparison of the outcomes of the Star lethal x 
Star, lethal crosses with those of the Star, Normal & Star, Normal crosses 
(table 1). The frequency of the wild-type class in the progeny of the 
former is lower (41.47 per cent) than in the latter (46.32 per cent). The 
wild-type flies in the progeny of Star/lethal X Star,/lethal crosses are 
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heterozygous for two (non-allelic) lethals, while their Star sibs are hetero- 


zygous for only one lethal (not counting Star itself). 

The second-chromosome lethals used in our experiments seem, on the 
average, to reduce the viability of the heterozygotes. This inference can 
be tested further by examining the frequencies of the wild-type class 
separately in various series of experiments summarized in table 1. The 


lethal and the Normal chromosomes were derived from four different 
localities Quarai, Catuni, Emboaba and Caxias. Seven possible com 
parisons can thus be made of the frequencies of the wild-type class in 
crosses which were free of lethals, and which contained lethals. These 
comparisons are as follows: 


STAR X NORMAI STAR X STAR 
LuTHAL LETHAL-FPREI LETHAL LE THAL-FREEK 


Quarai 60.42 60.33 $3.87 14.35 
Catuni 55.38 56.16 410 26 45.29 
Emboaba 51.29 55.94 44.87 19 45 
Caxias ie 41010 W315 


In six out of the seven comparisons, the lethal-containing wild-type class 
is less numerous than the lethal-free class in the corresponding crosses. 

In the comparisons just considered, the percentages of wild-type flies 
were computed from the total yield of the flies obtained in all the crosses 
of a given type. It will, however, be shown below that the outcomes of 
the crosses involving different combinations of the chromosomes are sig 
nificantly heterogeneous, as are different cultures involving the same com 
bination of chromosomes. Mean frequencies of wild-type flies were there 
fore calculated, first averaging the frequencies for the different cultures of 
each combination of chromosomes, and then averaging the average fre 
quencies for the different combinations of chromosomes. Since the num- 
bers of cultures, and of flies within the cultures, obtained for different com 
binations of chromosomes were unequal, these mean frequencies differ 
somewhat from the frequencies weighted in accordance with the numbers 
of the flies counted. The mean frequencies are shown in the rightmost 
column in table |. Mr. Allan Birnbaum and Prof. Howard Levene, of the 
Department of Mathematical Statistics, Columbia University, kindly 
devised a method for applying the analysis of variance to the data here dis- 
cussed. Table 2 shows the differences between the mean frequencies of 
the wild-type flies in the lethal-bearing and in the lethal-free progenies of 
the various crosses. The table also shows the standard deviation of these 
differences, from which ¢ values are derived. All the differences are nega 
tive, i.e., the lethal-bearing wild-type classes appear to survive less often 
than the corresponding lethal-free classes. The probability that all the 
differences would have the same sign by chance is 1/64. Considered 
separately, only two of the seven ¢ values are significant; however, these 





Vou. 38, 1952 GENETICS: A. R. CORDEIRO 175 


are also the only two based on enough data to have adequate sensitivity.’ 
Combined ¢ values for all the Star X Normal and all the Star & Star 
crosses are significant, and the combined ¢ value for all seven crosses, 
— 3.96, is highly significant. 

Estimation of the Magnitude of the Average Deleterious Effect Produced 
by the Lethals in Ieterozygous Condition. The difference between lethal 
free and lethal-bearing classes of flies is not an adequate measure of the 
relative viability of these classes, since this difference depends on the fre 
quencies of the Star and wild-type flies, and therefore on the viability of 
the carriers of the Star chromosome To measure the deleterious effect of 


heterozygosis for the lethals, a ‘coefficient of inviability’’* is used. This 


coetticient 1s not identical with the selection coeflicient (.S); it reasures the 


PABLE 2 
DIFFERENCE BETWEEN THE MEAN FREQUENCIES (IN PER CENT) OF THE WILD-TYPE 
CLASS IN THE LETHAL-CONTAINING AND THE LETHAL-FREE CROSSES 


STANDARD 
CROSS AND LOCALITY DIPPFERENCE DEVIATION 
Star * Normal 
Quarai —3.76 
Catuni —1.80 
1 23 


Emboaba 
Combined Star X Normal 
Star X Star 

Quarai 

Catuni 

Emboaba 

Caxias 
Combined Star & Star 
Combined all crosses 

“ Significant. 

” Highly significant 


proportion of the lethal-bearing individuals which are eliminated due to 
differential mortality between the egg stage and the hatching of the adult 
insect, in relation to the lethal-free individuals which survive. 

The coefficients of inviability were calculated only for the crosses 
Star lethal & Normal’) Normal and Star’ Normal X& Normal’ Normal. 
The other crosses, which contain two lethals, raise the problem of the 


possible interactions of the deleterious effects of heterozygous lethals. 


The coefficient proved to be 14.8 per cent for Quarai, 6.1 per cent for Catuni 
and 0.8 per cent for Emboaba, giving a mean estimate of 7.2 per cent. 
Tleterogeneity of Effects of Different Chromosomes... The data presented 
above show that, on the average, individuals heterozygous for lethals have 
a viability lower than individuals which are free of lethals. It is, however, 
probable that different lethals, as well as different chromosomes free of 
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lethals, affect the viability to different extents. An analysis of variance 
has been carried out under the direction of Mr. Allan Birnbaum _ Its 


results are summarized in table 5 

As indicated above, different lethal-bearing and lethal-free (Normal) 
chromosomes have been utilized in the experiments. <A total of S9 dif 
ferent combinations of chromosomes have been made in the crosses of the 
type Star/lethal & Normal/Normal (table 1). Some of these combina- 
tions have, however, been replicated in several culture bottles; a total of 


rABLE 3 
ANALYSIS OF VARIANCE OF FREQUENCIES OF WILD-TypE FLIES OBSERVED IN DIFFERENT 
CROSSES 


SUMS OF DEGREES OF MEAN 
TYPE OF CROSS SOURCE SQUARES FREEDOM SQUARES 


Star/lethal < Between chromo 
Normal/ Normal somes 13.132 x6 152.7 
Within chromo 
somes 11,246 ¢ 35 


Total 


Star/Normal Between chromo 
Normal / Normal somes 
Within chromo 


somes 


Total 
Star/lethal x Between chromo 
Star/lethal somes 
Within chromo 
somes 3461 


Total 5307 
Star/Normal Between chromo 
Star/ Normal somes 1516 
Within chromo 


somes 362 
Total 1878 30 


109 cultures were raised for crosses of this type. The analysis of variance 
shows that the variance between different combinations of chromosomes is 
significantly greater than that between cultures of the same combination 
of chromosomes. The F value equals 4.345, which is significant at a level 
below 0.5 per cent (table 3). Similarly, for the crosses of the type Star, - 
Normal & Normal, Normal, the F value is 2.112, which is significant at 
about the 2.5 per cent level (table 3). The crosses of the type Star lethal X 
Star/lethal and Star;Normal X Star/Normal do not give significantly 
greater heterogeneity between combinations of chromosomes than_ be- 
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tween replicate cultures. This is probably due to the small number of 
replications made in this series of experiments. 

This result is important. It shows that the so-called “Normal” flies in 
a natural population are in reality a diversified collection of genotypes hav 
ing different reaction norms and different adaptive values. The chromo- 
somes which carry recessive lethals, as well as chromosomes free of such 
lethals, contribute to the adaptive polymorphism of the normal species 
population. There is an indication in the data that different lethal-bear- 
ing chromosomes are more diversified with respect to their effects on 
heterozygotes than are the lethal-free chromosomes, but this point requires 
a more rigorous prool than can be adduced on the basis of the data so far 
available. There are also indications that while some lethal-bearing 
chromosomes are distinctly deleterious to heterozygotes, others are not 
deleterious, or even heterotic. It should also be noted that there exists a 
significant heterogeneity of the viability effects between the chromosomes 


derived from different localities (table 1). 

Acknowledgments... The writer wishes to express his appreciation to Mr. 
Allan Birnbaum, Prof. Howard Levene and Mr. Timothy Prout for ad- 
vice in the statistical treatment of the data; to Prof. Th. Dobzhansky for 
his help in the preparation of the manuscript; to Miss F. Wolff and Mrs. 
Malvina Cordeiro for assistance in the execution of the experiments. A 
part of the work has been done during the tenure of a Fellowship provided 


by the Rockefeller Foundation. 

Summary. The effects in heterozygotes of lethals found in second 
chromosomes of natural populations of Drosophila willistoni have been 
studied. On the average, the viability of lethal heterozygotes 1s_ sig 
nificantly lower than that of lethal-free individuals. The coefficient of 
inviability is of the order of 7 per cent. Both lethal-bearing and lethal 
free chromosomes vary significantly with respect to their viability effects 
in combinations with other chromosomes. 


1 Dobzhansky, Th., Genetics and the Origin of Species, 3rd ed., New York (1951). 

2 Berg, R. L., C. R. (Doklady) Acad. Sci. URSS, 3€, 212-218 (1942), 

3 Stern, C., and Novitski, E., Science, 108, 538-539 (1948) 

' Dubinin, N. P., Genetics, 31, 21-38 (1946). 

5 Wright, S., Dobzhansky, Th., and Hovanitz, W., /bid., 27, 363-394 (1942 

® Muller, H. J., Am. J. Human Genetics, 2, 111-176 (1950) 

7 Pavan, C., Cordeiro, A. R., Dobzhansky, N., Dobzhansky, Th., Malogolowkin, C., 
Spassky, B., and Wedel, M., Genetics, 36, 18-30 (1951). 

* Moree, Ray, Am. Natur., 86, 45-48 (1952) 

* Footnote by Howard Levene. A preliminary Chi-square test showed significant 
heterogeneity between replicate cultures. Because of this the usual Chi-square tests 
and the formula pqg/n for Dinomial variance are inapplicable and were never used. For 
all further work the percentage of wild flies per culture was treated as a normal variable 
with a variance to be found empirically. The analysis of variance in table 3 was an 
ordinary one-way analysis of variance with unequal numbers of replications, each 
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locality being treated separately and the sums of squares then combined. Since this 
analysis of variance showed significant differences between different combinations of 
chromosomes, for further work the average of the frequencies of wild-type flies for the 
different cultures of a single combination of chromosomes was considered as the basic 
variable. Then Star/Lethal X Normal/Normal at Emboaba, for example, might have 
a smaller frequency of wild type than Star/Normal X Normal/Normal at Emboaba 
simply because the chromosome combinations of the first type chosen by chance hap- 
pened to give a smaller frequency of wild-type. We can say the difference is due to the 
presence of the lethal chromosome only if it is significantly larger than the variability 
between chromosome combinations. This can be tested by an ordinary two sample 
i test asin table 2. Because of small numbers in many of the fourteen groups an esti- 
mate of the error variance based on all the data was used. Then ¢ is approximately 
normally distributed with unit variance, and the ¢ values can be combined by adding 
them (taking account of signs) and dividing by the square root of the number of t's 
added. This combined value will again be approximately normal with unit variance. 
Returning to the example of Emboaba, Star/lethal * Normal/Normal is based on 
twenty chromosome combinations and has adequate accuracy. Unfortunately, Star / 
Normal & Normal/Normal is based on only two combinations and this is the limiting 
factor in the sensitivity of the comparison. Most of the labor in obtaining the twenty 
combinations is wasted; in fact, three combinations of each type would have given a 
more accurate comparison. A more delicate and complicated test, based on the analysis 
of variance, gave substantially the same results. Mr. Birnbaum and the author of this 
footnote hope to discuss the statistical aspects more fully elsewhere. 


SELECTION AGAINST HETEROZYGOTES FOR) AUTOSOMAL 
LETHALS IN) NATURAL POPULATIONS OF DROSOPHILA 
WILLISTONI 


By Timotuy Prout 
CoLuMBIA UNIVERSI?Y, NEW YorRK 
Communicated by Th, Dobzhansky, March 31, 1952 


In an accompanying article, Cordeiro shows that the second chromosomes 
of Drosophila willistont which are lethal when homozygous, produce a 
significant lowering of the viability of heterozygous flies as well. The 
method used in the experiments of Cordeiro is comparison of the ratios of 
flies heterozygous for lethals and of flies which are free of such lethals. 
Cordeiro estimates that, under the environmental conditions of his experi- 
ments, the average lowering of the viability produced by heterozygosis for 
a second chromosome lethal is approximately 7°; (his coefficient of in- 
viability). The present note reports an attempt to compare the average 


adaptive values of individuals heterozygous for second chromosome 
lethals and of individuals free of such lethals in natural populations of D. 


willistont. 
The genetic structure of natural populations of D. willistoni has been 
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studied by Pavan, ef al.,| who have published data on the frequencies of 


second chromosomes which are lethal or semilethal to homozygotes. 
Dobzhansky, Spassky and Spassky* have obtained estimates of the muta- 
tion rates which produce these lethals and semilethals. Ives* and Wallace* 
have computed the probable number of loci which mutate to lethals and 
semilethals in the second chromosome of D melanogaster. The second 
chromosomes of D. willistont and D. melanogaster are homologous and 
probably contain the same loci.’ Professor C. Pavan* has kindly com- 
municated to us his unpublished data on the frequency of allelism of second 
chromosome lethals isolated from various Brazilian populations of D. wil- 
listont. 

The method used for the estimation of the adaptive value of hetero- 
zygotes for lethals in natural populations involves the comparison of the 
rates of the origin of such lethals by mutation and of their elimination in 
homozygotes. These rates should be equal if all lethals in natural popula- 
tions were completely recessive. Wright” * has given the necessary equa- 
tions which describe the relationship just stated (Wright's notation will be 
followed throughout). For the case of complete recessivity the balance of 
gains and losses at equilibrium is as follows: 


PO*X(1 — Q) 
1 — PO 


V — PQ 
where: 


Q = the equilibrium frequency of lethal chromosomes, 

V = the mutation rate producing lethal chromosomes, 

P = frequency of allelism of lethal chromosomes, 1.e., the proportion of 
zygotes bearing two allelic lethal chromosomes among the 
zygotes which carry two lethal chromosomes. 


The left side of equation (1) represents the rate of gain by mutation and the 
right side the rate of loss through homozygosis. Q, | and P may be esti- 
mated experimentally. It should be noted, however, that the experi- 
mentally determined value of P will be an underestimate if an appreciable 
amount of inbreeding occurs normally in the population. If the mutation 
rate (I’) is higher than the rate of elimination (PQ"), as is the case in all 
Drosophila populations studied so far, then the amount of selection against 
lethal heterozygotes required to account for this difference should be as 
follows: 


V — POU + V) 
QU. + 2V) — 2PQ7 + ‘V) 
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where: 


S = the mean selection coefficient of heterozygotes for lethal chromo- 
somes. 


Equation (2) still assumes random breeding. Wright’ * has considered 
the more general case where some elimination of lethals will be due to in- 
breeding as well as to semidominant deleterious effects of the lethals. 
There is evidence, however, that the populations of D. wiillistoni are 
panmictic,’ and that their genetically effective size may be taken as in- 
finite (unpublished data of Prof. C. Pavan). 

In such populations the frequency of the allelism is approximately equal 
to the frequency of allelism of independently arising lethal chromosomes 
(P.. according to Wright's notation). For the second chromsome of 
D. melanogaster,* * this value is P., = 0.0026 (+0.00057). Assuming the 
second chromosomes of D, melanogaster and D. willistoni carry the same 


TABLE 1 


FREQUENCIES OF LETHAL CHROMOSOMES (QV) AND SELECTION COEFFICIENTS OF HETER 
OZYGOTES FOR LETHALS (.S) IN VARIOUS POPULATIONS OF Drosophila willistoni 
POPULATION oO S 

Belem, Para 339 O25 
Rio Branco and Rio Negro 541 OLS 
Arce 13 OW 
Goyaz 417 020 
Catuni, Bahaa 143 O18 
Pirassununga, Sao Paulo 1253 O19 
Mogi, Sao Paulo 386 O21 
Vila Atlantica, Sao Paulo 392 O21 
Iguassu, Parana 375 022 

Mean 412 020 


equation (2), together with the estimate of the mutation rate, 17 = 0.0088, 
obtained by Dobzhansky, Spassky and Spassky.’* 

In table | there are presented the values of S obtained for nine natural 
populations of D. wil/istont in which Pavan, ef al.'! determined the frequen 
cies of lethal and semilethal chromosomes (semilethals are assumed to be 
effectively lethal under natural conditions). Diiferent values of P were 
tried out by substitution in equation (2) in order to determine how great 


are the differences in the estimates of S which depend on experimental 
errors in the determination of P. If values of P twice as high or half as 
high as 0.0026 are used, the resulting estimates of S are not appreciably 
changed. 

Considering the sources of error contributing to these estimates of S, a 
discussion of the variations of S between populations would be meaning- 
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less. All the examined populations of D. willistont may be roughly char- 
acterized by the value obtained using the mean Q, namely S = 0.020 
This value is of the same order of magnitude as the estimates made for 
natural populations of D. pseudosbcura,’* 0.014, 0.02, 0.007 and of D. 
melanogaster,*:* 0.054, 0.028, 0.016, O.OLL. It therefore seems justified to 
state that lethals in natural populations of Drosophila are on the average 
slightly deleterious in the heterozygous condition. 

The above estimates for natural populations vary from 0.007 to 0.05 for 
different species and different populations. The experiments of Cordeiro 
reported in the preceding article show, however, that under laboratory 
conditions the viability of heterozygotes for lethals in D. willistoni is 
lowered in some localities by as much as 14% (his coefficient of inviability), 
although in other localities it is only lowered by 0.8%. These estimates are 
not exactly comparable to the estimates of S calculated in this paper. 
In the first place, the experiment of Cordeiro measures the relative viabili- 
ties of different genotypes between the egg stage and the hatching ef the 
adults from pupae under conditions of severe crowding in laboratory 
bottles. These estimates of S measure the adaptive values under natural 
conditions where larval competition for food may not always be a limiting 
factor. On the other hand, the relative fecundities, longevities, sexual 
activities, etc., of the carriers of different genotypes influence the value of 
S, while they do not affect the viability estimates, such as those offered by 
Cordeiro. Finally it should be noted that the estimates of .S, as well as of 
the coefficient of inviability, describe merely the average situations. In 
reality, some lethals may be more deleterious 1n heterozygotes than others, 
as actually demonstrated by the data of Cordeiro. 

Acknowledgments. The writer wishes to express his appreciation for the 
valuable advice and suggestions of Prof. Th. Dobzhansky. This article 
was written during the tenure of a U. S. Atomic Energy Commission fel- 
lowship. 
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THI ACCUMULATION OF a-PHENYLGLYCINE BY MUTANTS OF 
NEUROSPORA CRASSA STIMULATED BY PHENYLALANINE 
AND TYROSINE* 

By C. H. Happoxt 


GENETICS LABORATORY, DEPARTMENT OF ZOOLOGY, THE UNIVERSITY OF TEXAS, AUSTIN, 
Texas} 


Communicated by J. T. Patterson, April 22, 1952 


During the course of experiments testing the effects of various muta- 
genic agents, two mutants of Neurospora crassa were obtained which grew 
slightly on the minimal medium but required phenylalanine or tyrosine 
for optimal growth. Genetic tests made it evident that these strains 
differed from the wild type by a single gene. It is the purpose of this paper 
to describe these mutants, one of which accumulates large amounts of 
a-phenylglycine; and to present evidence in favor of the hypothesis that 

TABLE 1 
GROWTH OF THE MUTANTS WITH PHENYLALANINE OR TYROSINE 


Recorded mg. dry weight (average of 2) of mycelium produced in 72 hours 
SUBSTANCE ADDED UT 4 UT 349 Em 5256 
Nothing 2.9 4.6 66.2 
Tyrosine concentration, molar: 

0. O00004 11. > 66.0 

0. 0000207 63.2 

0. 000040 19.9 68.2 

0. 000200 52 ' 64.2 
Phenylalanine concentratiop, molar: 

0 QO0004 5 64. 

0 000020 j : 69.2 

0 000200 é 65. 
Tyrosine and phenylalanine equimolar 

concentrations, molar: 

0 QO0004 12.0 28 | 

0. 000020 17.0 40 

0. 000040 22.3 46 

0, 000400 25.1 59 

® The value of 0.000020 molar, frequently used in other experiments, corresponds to 
an addition of 0.091 mg. of tyrosine and 0.083 mg. of phenylalanine per 25 ml. of medium. 


this compound may be utilized as a precursor of aromatic compounds 
under certain conditions. a-phenylglycine, a lower analog of phenyl- 
alanine, has heretofore not been considered a normal constituent of living 


matter. 

Evidence has been accumulating that indicates a common precursor for 
phenylalanine, tyrosine and other aromatic compounds. Haskins and 
Mitchell! have described a mutant of Neurospora that may utilize cinnamic 
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acid as an apparent precursor for phenylalanine and tyrosine or a common 
intermediate. Others have postulated that quinic acid or shikimic acid 
may act as a precursor for the synthesis of aromatic compounds in £s- 


cherichia colt and some Neurospora mutants.** There are also some indi- 
cations that the aromatic amino acid requirements of Leuconostoc mesen- 
teroides may be relieved by phloroglucinol.‘ 

Experimental.—The mutant strains, UT 39 and UT 4, were isolated as 
ascospores from crosses using ultraviolet-irradiated macroconidia (Em 
5297a) and untreated wild type (Em 5256A) protoperithecia by the method 
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DA YS 
FIGURE 1 


Growth of UT 39 on minimal and 2 * 10°° AM tyrosine and phenylalanine 
as compared to Em 5256A on minimal 


of Lein, Mitchell, and Houlahan.® UT 39 requires phenylalanine or tyro 
sine and UT 4 requires tyrosine for optimal growth. Growth was compared 
to Emerson Wild Type 5256A. 

Genetic investigations by Haddox® showed that the mutants UT 39 
and UT 4 are due to single gene differences. ‘The genes are not allelic, 
UT 39 being in linkage group D of Houlahan, ef a/.,? and UT 4 being in 
another, undetermined, linkage group. 

The minimal medium used in these studies was that of Beadle and Tatum. 
The growth was measured in terms of the dry weight of mycelium produced 


& 
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in 72 hours at 25°C. on 25 ml. of medium with certain exceptions as noted 

Growth Characteristics of the Mutants: UT 39 and UT 4 produce only a 
few milligrams of mycelium at the end of the third day of growth on mini- 
mal medium. Thereafter their growth rates increase greatly and by the 
eighth day they produce almost as much mycelium as wild type (Figs. | 
and 2). The addition of tyrosine or phenylalanine to the minimal medium 
stimulates the growth of UT 39, but not during the first two days (table 1, 
Fig. 1). When tyrosine and phenylalanine are present together in equi- 
molar concentrations, however, the growth approximates that of wild 


WILD TYPE MINIMAL 

UT4A — MINIMAL 

UT4A —MINIMAL + TYROSINE 
UT4A -~MINIMAL + @-AL ANINE 


GROWTH CURVE 
UT4A 


FIGURE 2 


6 7 











DAYS 
FIGURE 2 


Growth of UT 4 on minimal and 2 * 10°° Mf tyrosine and phenylalanine 
as compared to Em 5256A on minimal, 


type after three days (table 1). UT 4 is stimulated only by phenylalanine, 
and the presence of both amino acids together has no significant effect 
on its growth rate. Neither mutant is a pH or temperature sensitive 
mutant. 

All known amino acids, vitamins and purines and pyrimidines as well 
as a number of organic acids and other compounds known to have growth- 
promoting properties were tested to determine their effect on the mutants. 
In addition, iodoacetate, KCN, sulfasuxadine and taurine were tested for 
possible inhibition. Some of these showed an inhibitory or stimulatory 
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effect particularly in the presence of tyrosine or phenylalanine. These 
effects are summarized in table 2 in which stimulation or inhibition is 
given as percentage of growth of the mutant calculated on the basis of wild 
type growth on minimal medium as 100°). It will be noted that none of 
the compounds tested stimulated mutant growth as much as tyrosine or 
phenylalanine. Two compounds, leucine and 3-hydroxyanthranilie acid 
showed an unexpected inhibition. Tryptophan is known to be slightly 


inhibitory to wild type, and therefore its effect on the mutants is not sur- 
prising. 

It was found also that iodoacetate inhibited both the mutants and wild 
type at concentrations of 0.001 \/. The inhibition is partially reversed 


TABLE 2 
GROWTH PROMOTING OR INHIBITING SUBSTANCES 


Values listed are dry weights of mutant and wild type mycelia, expressed as percentage 
of the average wild type mycelial dry weight when grown on minimal medium for 72 
hours at 25°C. 

UT 4 UT 39 
PLUS PLUS PLUS PLUS Sas 
TYRO PHENYI T¥YRO PHENYL 5256, 
SUBSTANCE ADDED ALONE SINK@ ALANINE ALONE SINE ALANINE ALONE 
Nothing 39.1 ‘4 5.1 33.$ 3 100.0 
3-Hydroxyanthranilic acid 
(40 yg.) : oa hf 114.0 
Leucine (1.0 mg.) > s ‘ ae § ie OS 
Tryptophan (1.0 mg.) 24 a! i 23 2 ; 75.8 
a-Aminobutyric acid (0.5 
mg.) < 658 2.% i} ee f d 99 
Citric acid (2.5 mg.) 10 78 5.6 6 { 108 
Coenzyme I (0.1 mg.) 18 67.3 3.6 42.5 47 ¢ - 102 
Fumaric acid (2.5 mg.) 10.: 64.5 , 14.; ii 100 
a-Ketoglutaric acid (0.5 
mg.) 14.3 63.5 16 32.0 5S 100 
Malic acid (2.5 mg.) 2.6 70.0 & 17 4 ; 78 115.0 
Succinic acid (2.5 mg. ) 17.0 84.0 14.3 26.7 7 85 112.0 
d-Phenylalanine (O.0l mg.) 3.0 31.0 a 34.0 78 7TH0 100.0 
“ The amounts of tyrosine or phenylalanine added are those considered to be necessary 
for approximately half maximum growth, or 2 * 10 ° molar 


by tyrosine, phenylalanine or succinic acid (table 3 ) in all three strains, 
but whereas other amino acids are equally effective for relief of inhibition 
in wild type, only phenylalanine and tyrosine are effective for the mutants. 

Accumulation of an Active Growth Factor: Bacterial assays of the culture 
media in which UT 39, UT 4 and Em 5256 had grown for a period of three 
days indicated an apparent accumulation of phenylalanine and tyrosine 
by the mutants but very little by the wild type. Using Leuconostoc 
mesenteriodes P 60,and Lactobacillus arabinosus USDA 17-5 afterthe method 
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of Hac, Snell, and Williams,’ it was found that UT 89 growing on minimal 
medium accumulated after 72 hours 22.5 wg. per mg. dry weight of myce- 
lium of a tyrosine active compound and 7.6 yg. of a phenylalanine active 
compound. Comparable figures for UT 4 were 2.5 ug. for tyrosine and 
13.8 wg. for phenylalanine. Wild type produced | yg. or less in both assays 
under the same conditions. 

In order to supplement the data obtained with the bacterial assays, the 
culture medium was chromatographed on paper. No tyrosine or phenyl- 
alanine spots could be developed with ninhydrin, but a spot did appear which 
in area was roughly proportional to the apparent phenylalanine and tyro- 
sine concentrations as determined by the biological assays. The conflict- 
ing results of the bacterial assays and the paper chromatograms suggested 
the accumulation of some factor which satisfied the requirements of the 
Neurospora mutants as well as the bacteria for phenylalanine and tyrosine. 

Cross feeding experiments showed that both UT 4 and UT 39 produce a 
compound or compounds in the medium after 3 days growth on minimal 
which stimulates their growth. UT 39 produces a highly active medium 
for both itself and UT 4, while UT 4 produces a relatively low activity 
compared to UT 39. Wild type medium was found to be completely inac- 
tive for either mutant, and the active factor produced by the mutants had 
no effect on wild type growth. 

Isolation and Identification of the Active Factor: The paper chromato- 
grams of the media from UT 39, UT 4 and Em 5256 showed a ninhydrin 
positive spot which had a different Rf value from the known amino acids. 
This substance was accumulated in large amounts in UT 39 medium only, 
and therefore this medium was used for its isolation. 

UT 39 was grown for eight days in minimal medium, the medium filtered 
to free it of mycelium and concentrated twentyfold in vacuo at ca. 60°C. 
Ethanol was added to 90% by volume to precipitate salts and sugars, 
and the active unknown concentrated on filter paper by using the chroma- 
topile technique of Mitchell and Haskins" with a water-saturated phenol 
solution as a solvent. ‘The location of the active factor in the filter paper 
pile was determined by running supplementary paper chromatograms and 
using the Neurospora mutants and Leuconostoc mesenteroides as assay 
organisms. The active factor was removed from the paper with ethanol 
and the alcohol and residual phenol removed by distilling 7m vacuo and re- 
placed with water. The pH of the water solution was adjusted to 7.00 with 
NaOH and an active crystalline precipitate was formed. This was washed 
with diethyl ether, 9007 ethanol and water, redissolved in NaOH and the 
solution filtered. The solution was again returned to neutrality, filtered 
and the crystals washed with 95°, ethanol and water. From 5500 ml 
of original culture medium in which 23.32 g. of UT 89 mycelium had ac- 
cumulated, 370 mg. of these purified white crystals were obtained. 
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TABLE 3 
IODOACETATE EFFECTS 
Recorded mg. dry weight of mycelia produced in 72 hours 
UT4 UTSY 
PLUS PLUS 
lobo 1ODO 
SUBSTANCE ADDED ALONE ACETATE@ ALONE ACETATE ALONE ACETATE 
Nothing §.2 0.0 6.2 0.0 64.0 7 
Tyrosine (0.001 molar) 24.0 10.0 31.0 5: 58.0 20.2 
Tyrosine (0.010 molar) 26.0 14.0 30.0 14.0 54.0 23 
Phenylalanine (0.001 molar) 1.3 5.0 28.0 18.0 58 22 
Phenylalanine (0.010 molar ) 6 32 24.0 56.0 24 
Succinie acid (0.001 molar ) 18.: 8.: 14.8 10.0 60.0 14 
Succinie acid (0.010 molar) 17.6 10.6 18 12.0 58.0 10 
Succinie acid (0.100 molar ) 20 18 22.5 17.8 60.0 48 


“The concentration of the iodoacetate was 0.001 molar. The iodoacetate 
sterilized by passing through a VF sintered glass filter 


TABLE 4 
Rf VALUES OBTAINED FROM UNKNOWN, d/-a@-PHENYLGLYCINE, AND CONTROLS 


Average of three determinations 


SOLVENT 
BUTANOL BUTANOL PHENOL WATER 
ACETIC ETHANOL ACKTIC SATURATED 
SUBSTANCE METHANOL ETHANOL ACID AMMONIA ACID PHENOL 


Aspartic acid 0.24 0.05 26 03 0.25 0.05 
Phenylalanine 0.73 0.50 69 0.44 ‘ 0.76 
Tyrosine 40 0.19 16 0.18 35 0.50 
Tryptophane 60 0.16 54 0.40 7 0.67 
a-Phenylglycine 92 0.73 71 0.34 ae 0.19 
Unknown 92 0.75 71 0.33 76 0.19 





TABLE 5 
GROWTH STIMULATION BY UNKNOWN AND dl-a@-PHENYLGLYCINE 


Recorded mg. dry weight of the mycelia produced in three days. The Leuconostoc 
readings are visual estimates of proportionate growth volume 





Leuconostoc 
mesenteroides 


~ 


SUBSTANCE ADDED 
Tyrosine (0.02 mg.) 
Tyrosine (0.10 mg.) 
Phenylalanine (0.02 mg.) 


KS # b> 


Phenylalanine (0.10 mg.) 
Unknown (0.02 mg.) 

Unknown (0.10 mg. ) 14 
dl-Phenylglycine (0.02 mg. ) 5.2 
dl-Phenylglycine (0.10 mg.) 17.4 
Nothing ‘ 1.0 


te & be te 


25 
2 


“ Basal medium supplemented with phenylalanine 
* Basal medium supplemented with tyrosine. 
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The crystals decomposed between 230 and 260°C., showing no melting 
point. The compound was found to react with pyruvic acid in the Strecker 
reaction'! with the production of alanine (identified by paper chromatog- 
raphy) and benzaldehyde. The benzaldehyde was identified by precipi- 
tation with phenylhydrazine to give the reddish white crystals of the phenyl- 


hydrazone of benzaldehyde (m.-p. 154-157°C.). 

From this it was indicated that a-phenylglycine might be the substance 
accumulated. Synthetic d/-a-phenylglycine was prepared after the method 
of Steiger.'” Paper chromatograms were run on both the unknown and the 
synthetic d/-a-phenylglycine. The spots are identical (table 4). Growth 
experiments with the unknown and d/-a-phenylglycine indicated that the 
effects on growth are identical (table 5). 

Discussion..-The above data show that the phenylalanine and tyrosine 
requirements of two mutants of Neurospora crassa and a strain of Leu- 
conostoc mesenteroides may be replaced by a-phenylglycine. One of the 
mutants, UT 39, accumulates considerable amounts of this substance 
Whether this should be considered a normal precursor, is open to question. 
There are at least two possibilities. First, it may be a normal precursor 
produced in the cell and normally immediately utilized. A second pos- 
sibility is that it is a product of “shunt metabolism.’’ Thus the normal 
synthesis of phenylalanine and tyrosine is blocked in UT 39 allowing the 
accumulation of an unstable precursor which is converted to a-phenylgly- 
cine. 

The sudden increase in growth of the mutants after a period of days is 
best explained by postulating that an enzyme system, already present but 
inefficient in the utilization of a-phenylglycine, is slowly adapted to more 
efficient usage of available a-phenylglycine. Once it has been adapted, 
it may perpetuate itself only in the presence of the a-phenylglycine. 

The accelerated growth of the mutants after a period of time, in the 
absence of tyrosine or phenylalanine, resembles certain of the findings of 
Simmonds!® with Escherichia colt and Lactobacillus arabinosus. It is 
possible that the stimulation that she found with the bacteria may be due 
to the elaboration of a substance identical with that produced by the mutants 
used in these studies. 

The stimulation of the mutants by intermediates of the Kreb’s tricar- 
boxylic acid cycle, the stimulation by Coenzyme I, and the reversal of 
iodoacetate poisoning by tyrosine and phenylalanine may or may not 
indicate that these compounds are involved directly in the formation and 
utilization of a-phenylglycine. A number of hypotheses may be made to 
explain their stimulatory effect, but it is believed that the data are insuffi- 
cient at the present time to present a workable hypothesis. 

The reason for stimulation of the mutants by a-aminobutyric acid is 
unknown. A possible explanation is that a-aminobutyric acid functions 
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directly in the formation of one of the aromatic amino acids or in the 
conversion of a-phenylglycine to tyrosine or phenylalanine. A more in- 
direct function, however, is equally possible on the basis of the data pre- 


sented. 

It should be noted in conclusion that UT 39 when combined with the 
pantothenicless mutant 5531 to form the double mutant relieves the re- 
quirement for pantothenic acid.'* The reason for the suppression of the 
pantothenate requirement by a mutant gene causing the accumulation of 


a-phenylglycine is not obvious. 

The author would like to acknowledge the guidance given by Dr. R. P. 
Wagner in these studies. 

Summary.—Mutants of Neurospora crassa requiring tyrosine or phenyl- 
alanine for optimal growth were studied. One of these mutants produces 
an excessive amount of a-phenylglycine which under certain conditions 
may replace the aromatic amino acid requirements of the mutants and 
Leuconostoc mesenteroides. The significance of this compound, which is 
here for the first time proved to occur naturally, is discussed. 


* Taken, in part, from a thesis presented to the Graduate School of the University of 
Texas in partial fulfillment of the requirements for the degree of Doctor of Philosophy. 

+ Present address: Department of Surgery, Louisiana State University School of 
Medicine, New Orleans, La. 

t A portion of this work was conducted at the Department of Surgery, Louisiana State 
University School of Medicine, New Orleans, La. 
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THE EFFECT OF OXYGEN CONCENTRATION ON 
X-RAY-INDUCED CHROMOSOME BREAKAGE IN MAIZE* 


By DREW SCHWARTZ 
BroLoGy Diviston, OAK RIDGE NATIONAL LABORATORY, OAK RIDGE, TENNESSEE 
Communicated by M. M. Rhoades, March 24, 1952 


The effect of reduced oxygen tension in decreasing the radiosensitivity 
of biological material to x-rays has been demonstrated in many organisms. 
Experiments with onion,' Tradescantia,?~* Drosophila® and barley® have 
established that the production of x-ray-induced chromosomal aberrations 
is influenced by the amount of oxygen present in the cell at the time of ir- 
radiation. There are two ways by which the frequency of aberrations can 
be reduced, (1) by a reduction in the number of effective ionizations and (2) 
by an increase in the rejoining power of broken ends. The results of earlier 
investigations have been interpreted as favoring the first hypothesis. It 
was postulated that x-rays produce chromosome breaks in two ways, by 
ionization of the molecules of the chromosomes themselves and indirectly 
by ionization in the vicinity of the chromosome. In the latter case highly 
active intermediate products such as free radicals or peroxides are liberated 
which are capable of diffusing to the chromosome and causing breakage. 
Chromosomes irradiated in the presence of oxygen are thought to be broken 
by both mechanisms, while in reduced oxygen tension, the indirect effect 
is minimized. 

The experiments described here were conducted to determine if oxygen 
affects the radiosensitivity of maize chromosomes. Evidence will be pre- 
sented supporting the hypothesis that the effect of oxygen tension is on the 
recovery process rather than on the initial breakage mechanism. 

Experimental Methods. Genes a (aleurone color) and sh (shrunken endo- 
sperm) are situated on the long arm of chromosome 3. They are very 
closely linked giving 0.27°( crossing over.’ Pollen carrying both dominant 
genes was x-rayed at 1000 r and used to pollinate tester plants homozygous 
recessive for both genes but carrying all the other factors required for the 
expression of color in the aleurone, A,, Cand R. The exposures were made 
with a Collidge self-rectifying tube with a tungsten target, operated at 
250 kv. and 15 ma. (3 mm. of Al, inherent filtration). 

The pollen to be treated was collected from a number of plants, mixed 
thoroughly, and divided into two equal parts. One-half the pollen was ir- 
radiated in air and the other half in a 100°% nitrogen atmosphere. The 


exposures were made in a lucite chamber. In the case of the nitrogen 


treatment, the chamber was made airtight and the pollen preflushed with 
the gas for 3 minutes. A steady flow of nitrogen was passed through the 
chamber during the 4'/) minutes of exposure. 
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Two classes of mutant kernels were observed, whole seed losses where the 
entire kernel showed only the recessive phenotype, and mosaics which were 
sectorial for areas of both dominant and recessive tissue. Approximately 
3% of the kernels in the former class showed loss of only one of the two 
genes. All the kernels in the second class were mosaic for both genes. 

Since closely linked genes were used in this study it is possible to distin- 
guish the genetic effects due to chromosomal aberrations from those due 
to gene mutation. Losses of both A and Sh are due to chromosome break- 
age. Stadler and Roman* have shown that x-rays produce very few, if any, 
true gene mutations. Thus the chances that both genes A and Sh mutate 
simultaneously in the same pollen grain are practically nil. The single gene 
mutants are due to gene mutations or an internal deletion where one of the 
breaks occurred in the region between the two genes. 

Results and Discussion._-A comparison of irradiation effects produced in 
nitrogen and in oxygen shows a decrease in the number of mutant kernels 
by a factor of 3.2 (table 1). This over-all protection value compares well 
with that obtained with Tradescantia and Drosophila. The data also indi- 
cate that there is a significant difference between the relative protection of 





TABLE 1 
FREQUENCY OF MUTANT KERNELS FOLLOWING IRRADIATION IN AIR AND NITROGEN 


NO. OF WHOLE LOSSES MOSAICS TOTAL, 
SEEDS NO % NO \ % 













Air 6158 159 2.58 40 0.65 3.23 
Ne 4505 43 0.95 3 0.06 1.02 
Air/N» on 2.7 9.9 3.2 





low oxygen concentration in the two classes of mutant kernels. The 
whole losses show only 2.7-fold protection while the mosaics show 9.9-fold 
protection. This difference is significant at the 2% level (x* = 5.7) as 
determined by the likelihood ratio test and may offer a clue to the mecha- 
nism of oxygen protection. If reduced oxygen tension acts only to reduce the 
number of chromosome breaks, as has been postulated, there should be no 
difference in the direction found in the ratio of whole losses to mosaics 
when the exposures are made in oxygen or nitrogen. That a difference of 
this kind does exist suggests that it is the mechanism of chromosome resti- 
tution following breakage which is sensitive to the oxygen tension in the 
cell at the time of irradiation. 

This will be clarified by an analysis of the types of chromosomal aberra- 
tions which give rise to the two classes of mutant kernels. Let us consider 
first breaks limited to the long arm of chromosome 3. The mosaics arise fol- 
lowing a break distal to the A-Sh/ region. After chromosomal duplication 
the two broken ends of the sister strands fuse in endosperm tissue. Due 
to the chromatid type of breakage-fusion-bridge cycle a mosaic phenotype 
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results in the endosperm.’ The whole-loss kernels can arise from two types 
of aberrations, (1) chromosome breakage proximal to the A-Sh region re- 
sulting in a terminal deletion which includes the A and Sh genes and (2) 
chromosome breakage on both sides of the A-Sh region followed by recom- 


bination of broken ends in such a way as to result in an interstitial deletion. - 


It is apparent that the mosaic class results from chromosome breakage 
without reunion (with the exception of a few rings and dicentrics), while a 
portion of the whole-loss class (i.e., interstitial deletions) results from break- 
age with reunion of broken ends. If irradiation in nitrogen increases the 
ability of broken chromosome ends to rejoin, the number of terminal dele- 
tions, and hence dicentric chromatids, would be reduced but the interstitial 
deletions would not be reduced, since they represent a class where chromo- 
somal rejoining had, in fact, occurred, This would result in a more marked 
protection of low oxygen tension against the mosaic class as compared to the 
whole losses. 

Mosaic kernels may also arise from other aberrations such as centric 
rings and dicentrics produced by the fusion of the broken ends of two cen- 
tric chromosomes. These aberrations produce the chromosomal type of 
breakage-fusion-bridge cycle.'° In this study it was not possible to dis- 
tinguish between endosperm mosaics due to the chromatid and chromosome 
types of cycles. The frequency of occurrence of rings and dicentric chro- 
mosomes which contain the A-S/ region was determined in another set of 
experiments involving sporophyte tissue. The marker used in these stud- 
ies was the minute, cytologically undetectable a-x, deficiency found by 
Stadler and Roman which includes the A, Sh, and W (a chlorophyll gene) 
loci.® The chromatid cycle persists only in the gametophyte and endo- 
sperm while the chromosomal type continues in the sporophyte or plant 
tissue as well.'° Hence, plants mosaic for anthocyanin and chlorophyll 
pigmentation must be due to the chromosomal breakage-fusion-bridge 
cycle. It was found that a dose of 1000 r in air yielded 4 mosaic plants in 
a population of 5718. Assuming that the frequency of occurrence of ring 


and dicentric chromosomes is the same in the endosperm and the sporo- 
phyte 4 of the 40 endosperm mosaics obtained in the previous set of ex- 
periments should be due to these aberrations and 36 to terminal deletions. 
On the hypothesis that the oxygen effect is exerted on the recovery process, 


the kernels which are mosaic because of ring and dicentric chromosomes, 
would not be protected. Hence, the frequency of their occurrence should 
not be reduced, but perhaps increased, in the nitrogen treatment. Three 
such mosaics would be expected out of the 4505 kernels in the nitrogen- 
treated material. This would account for all the mosaics actually found. 
Thus, there appears to be an even more striking reduction in the occurrence 
of terminal deletions than the data actually show, 








ee 
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Breakage involving other chromosomes in addition to chromosome 3 does 
not affect this picture. Let us consider a case of breakage involving the dis- 
tal region of chromosome 3 and another chromosome. Reunion of the two 
centric fragments to form a dicentric chromosome is infrequent and has al- 
ready been considered. Rejoining of the chromosome 3 centric fragment 
with either acentric fragment results in a phenotypically normal kernel. 
Only those cases where the broken end of the No. 3 centric fragment re- 
mains free result in mosaics. Hence here, also, mosaics result from break 
age without reunion. Where the break in chromosome 3 occurred in the 
region proximal to the A-Sh loci only those cases where the acentric No. 3 
fragment failed to rejoin with either centric fragment give rise to whole 
losses. This class of whole losses is reduced in nitrogen treatment as are 
the mosaics. 

Further work is in progress utilizing a chromosome marked throughout 
its length, which will facilitate distinguishing between interstitial and ter- 
minal deletions. 

If oxygen is acting on the recovery process it must be an immediate ef- 
fect, induced at the time of irradiation, and which persists once it is ini- 
tiated. Giles and Riley* have shown the oxygen tension existing in the cell 
at the time of irradiation to be the important factor. Altering the oxygen 
tension immediately after irradiation does not alter the frequencies of 
chromosomal aberrations. Thus it must be postulated that there is a 
qualitative difference between broken ends of chromosomes, depending 
upon the amount of oxygen present at the time of breakage. Chromosomes 
broken in the absence of oxygen are more capable of rejoining than chromo- 
somes broken in the presence of oxygen. McClintock has presented evi- 
dence that the behavior of broken ends is ‘probably related to the method 
by which the chromosome becomes broken and to physiological conditions 
surrounding the broken end.’’® 

The data presented are not in agreement with the alternative hypothesis 
that fewer breaks are produced by irradiation in nitrogen. If this were the 
case the frequency of mosaics to whole losses should be higher in nitrogen 
than in oxygen. A reduction in the number of breaks would lessen the 
probability of the occurrence of two breaks in a single chromosome arm, 
thereby decreasing the frequency of interstitial deletions and hence whole 
losses. This is the reverse of what has actually been found. 

Summary.—-Experiments performed on the irradiation of pollen in air 
and 100% nitrogen atmospheres have demonstrated the effect of low oxy- 
gen concentration on the radiosensitivity of maize chromosomes. The dif- 
ferential protection of low oxygen tension against interstitial and terminal 
deletions suggests that the effect of oxygen is exerted on the recovery proc- 
ess rather than on the initial breakage mechanism. 





494 ASTRONOMY: W. J. LUYTEN Proc. N. A. S 


Acknowledgments.--The author wishes to acknowledge the assistance of 
Miss Rachel Clark. He is also indebted to Dr. A. W. Kimball for aid in the 
statistical analysis of the data. 

* Work performed under Contract No. W-7405-eng-26 for the Atomic Energy Com- 
mission 

' Thoday, J. M., and Reed, J., Nature, 160, 608 (1947) 

2 Giles, N. H., Jr., and Riley, H. P., these PROCEEDINGS, 35, 640 (1949) 

* Giles, N. H., Jr., and Riley, H. P., /bid., 36, 337 (1950). 

‘Giles, N. H., Jr., and Beatty, A. V., Science, 112, 643 (1950). 

5 Baker, W. K., and Sgourakis, E., these PROCEEDINGS, 36, 176 (1950). 

® Hayden, B., and Smith, L., Genetics, 34, 26 (1949). 

7 Mains, E. B., J. Heredity, 40, 21 (1949) 

* Stadler, L. J., and Roman, H., Genetics, 33, 273 (1948). 

* McClintock, B., /bid., 26, 234 (1941) 

” McClintock, B., these PROCEEDINGS, 28, 458 (1942) 


NEW WHITE DWARFS IN THE SOUTHERN HEMISPHERE: 
SECOND LIST 


By WiL_eM Jj. LuyTeNn 
UNIVERSITY OF MINNESOTA 
Communicated by Harlow Shapley, April 22, 1952 


Following up the color survey of proper motion stars in the Southern 
Hemisphere reported on previously,' results have now become available 


° 


for the zone with declination limits —55° and —60°. In the original 
survey for the detection of proper motions 48 pairs of plates taken with 
the Harvard 24-inch Bruce telescope were used. These were blinked by 
the writer and some 6500 stars were found to possess appreciable proper 
motion. These motions have all been published in part C of the General 
Catalogue of Proper Motions. Plates were taken for all 48 regions with 
the same telescope and on a yellow-sensitive emulsion; these have now 
been compared with the original blue plates and approximate colors for 


the 6500 proper motion stars determined. 

The vast majority of these stars are ordinary red dwarfs, of course, but 
some 250 stars were singled out as definitely bluer than expected. Narrow- 
ing down the limits still further, forty-three stars have been selected from 
among these 250 as holding out the greatest probability of being genuine 
white dwarfs. Data for these stars are given in table 1, the several columns 
of which require no further explanation. 


The forty-three stars have been arranged in four groups, as follows: 
(a) The first group contains the five stars which had previously been 
identified as white dwarfs from color plates taken at the Cordoba Observa- 














Dt 





VoL. 38, 1952 








BPM L 











16115 170-27 
16139B 170-14B 
18164 182-61 
26734 210-160 
27273 212-19 





16421 171-50 
16571 172-15 
17088 175-34 
18394 184-75 
18615 185-53 











18972 186-120 





19652 188-86 
19738 189-36 
} 21641 195-29 


22038 196-234 











200-115 





23550 


24047 202-179 
24723 203-131 
24866 204-118 
26536 209-44 


210-106 


214-57 









26744 
27891 








170-31 





16085 


16715 173-50 
19929 190-21 
20372 191-28 
21018 192-30 








21094 192-81 

22046 196-10 
| 22855 199-147 
23593 200-34 
24290 202-131 


24342 203-218 
24537 203-213 















204-16 





24919 


26691 210-114 
26998 211-90 
27561 213-30 








173-64 





16628 


16964 175-70 
} 180-12 
26183 208-31 
j 26658 210-25 











I 
20 1 
1 


21 


6 


9 34... 


50 


12 53.-+ 


13 27 


16 16. 
17 09. 
26. 


20 00, 


0 20 


2 00. 
10 13. 


59) 


11 29 


13 27. 
14 41. 
15 36. 


16 34 


OO .« 


20 14. 
40, 
21 46.: 


ro 


rh — 


§ 22 


19 27. 


20 10. 


xe 


1. 


oe. 


aad 


R.A. 1950 pec. 


— 55°42’ 


— 54 
—59 
— 58 
— 56 
— 58 
— 55 
— 56 
— 58 


iF 





— 59 


— 56 


—59 ¢ 


—57 : 


— 56 


—58 ¢ 


— 58 
—59 
— 56 


ad 
-— 


TABLE 


58 
11 
31 
03 
OS 
58 
34 
46 


02 


32 
43 
12 

09 


56 


09 
34 
53 


12 


18 
18 


04 
18 
49 


9 13 


16 


14.§ 


16 


15 
15 


14 


14 
13 
15 
14 


15 


15 
14 
16 
14 


14 


12 


13.§ 
13.¢ 


14 


15 
16 
16 


ASTRONOMY: W. J. LUYTEN 















aS e 


“NI te 


x 


uo & 























“ p Lc, 













0.585 214° 0"0 
0.340 242 +0.3 
0.301 184 +0.2 
0.511 144 —(.1 
0.450 114 +0.2 























0.122 99 +0.2 
0.155 100 +02 
0.118 76 0.0 
0.187 343 0.0 
0.074 336 0.0 











194 





298 +0. 








0.100 183 —(0.2 
0.188 140 0.0 
0.126 262 +0.3 
0.116 276 +-0.2 














131 219 +01 
0.243 209 —0.3 







0.205 207 0.0 
0.110 9 +0.2 
0.083 133 +0. 1 
















125 309 +0 .2 
200 96 0. 
















0.096 206 +().2 
0.095 193 +(),2- 
0.136 298 +0.2 
0. O85 274 +(0).2 
0.040 261 +01 























092 310 +0). 1 
0.128 246 +0.2 
0.101 156 +0.2 
0.065 211 +0.2 
0.063 178 +0.2 















O84 198 +O. 
0. O87 198 +03 




















110 216 +0). 
0.076 118 +01 
0.050 126 +0.3 
0.095 146 +03 



















0.051 141 






0.052 201 Less 
0.102 128 than 
0.031 206 | +0.4 






0.199 169} 


496 CHEMISTRY: J. E. LENNARD-JONES Proc. N. A. S. 


tory. Four of these were again recognized as white in the present survey 
but the fifth, L 170-14 B, is the faint companion of a close binary not 
resolved on the Harvard plates. 

(6) The second group comprises the seventeen stars for which the 
probability is very high that they are white dwarfs. 

(c) The third group contains sixteen further stars; probably at most 
one-third of these will ultimately prove to be white dwarfs. But while 
this will remove several stars from the list individually, others from among 
the 250 are expected to take their place and statistically speaking the stars 
in this group should be added to the number in group ()). 

(d) The fourth group contains those five stars which were invisible on 
the yellow plates. Past experience has shown that the chance of their 
being really white is very high but not much more can be said at the 
present time. 

It is a pleasure to acknowledge the financial support received from the 
Rumford Fund of the American Academy of Arts and Sciences, which 
made possible the taking of the yellow plates. The author is further 
indebted to the Harvard Observatory, and to its director, Dr. Harlow 
Shapley, for his cooperation in having these plates taken and his continued 
interest in and support of this program of research. 
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THE SPATIAL CORRELATION OF ELECTRONS IN ATOMS AND 
MOLECULES II* 
By J. E. LENNARD-JONES 
DEPARTMENT OF THEORETICAL CHEMISTRY, CAMBRIDGE, ENGLAND 
Communicated by R. S. Mulliken, March 31, 1952 


1. Introduction.-A discussion of the spatial correlation of electrons of 
the same spin in atoms and molecules has been given by the author in 
another paper presented to this Conference.' It was there shown that the 
exclusion principle is a potent factor in keeping electrons of the same spin 
apart. But to the approximation that the wave function of a many- 


electron system is represented by a single determinant there is no correla- 


tion between electrons of opposite spin. On physical grounds, however, 
we must expect some correlation between electrons, whether of the same 
or opposite spin, because of the electrostatic forces between them. This 
correlation is neglected when electrons of opposite spin are described by a 
simple product of two similar functions, each involving the coordinates of 
one electron only. This is done, for example, when the helium atom is de- 
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scribed approximately by a product of wave functions of the Hartree type, 
or when the electrons of the hydrogen molecule are described by a product 
of simple molecule orbital functions. In each case the distribution in space 
of one electron is assumed to be independent of the position of the other 
(except for some average of the electrostatic field). 

Attempts to improve this simple approximation to the wave function of 
a pair of electrons of opposite spin have been made by many authors by 
the method of ‘configuration interaction.’’ Essentially this method is an 
attempt to take into account the effect of electrostatic repulsion on spatial 
correlation, though its physical significance is not always realized nor are its 
implications discussed. Sometimes the correlation of electrons due to 
electrostatic repulsion is represented empirically by the introduction of the 
coordinate r,. into the wave function, as Hylleraas did for the helium atom 
and James and Coolidge*® for the hydrogen molecule. When this is done, 
the improvement in the calculation of the energy of the system is most 
marked, but the method has not proved to be popular because of the 
awkward form of the wave function and the difficulty of visualizing its 
properties. 

An attempt is made in this paper to introduce the effect of the electro- 
static field on spatial correlation in a more direct and logical way. The 
mutual repulsion of two electrons depends only on their distance apart but 
since the probability of the occurrence of the electrons in prescribed ele- 
ments of volume will, in general, be a function of six-coordinate space, it is 
desirable also to consider the electrostatic repulsion as a function of six 
dimensional space. The form of this function is well known when both 
electrons are described in terms of spherical polar coordinates referred to 
the same center. It is also known when the electrons are described by ellip- 
tic coordinates, as is convenient for two center problems. In each case the 
function is expressed as a series of terms, each of which has the same sym- 
metry in three-dimensional space but has a different symmetry in six- 


dimensional space. The wave function of the two-electron system has ac 


cordingly also been expanded in a series of terms with corresponding sym- 
metry properties. It is then found that the wave equation can be separated 
into a number of component parts, each of which has the same symmetry 
in six-dimensional space. This leads to a simplification of the general 
problem and to a closer understanding of the way the electrons affect each 
other in their spatial distribution. 

2. Helium and Similar Atomic Systems in the Ground State.—For a sys- 
tem of two electrons, each of which moves in the same central field V(r), 
the wave equation in atomic units 1s 


IT, + IT, of | ryo)W(1, 2) = Ey, 2), (2.1) 


where // as usual includes the kinetic energy operator and V(r). The ex- 
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pansion of 1/rj in spherical polar coordinates is 


1/rie = > UP, (cos 12), 


where 


U, = ri/r', (2.3) 


rT 


when 7; is smaller than r2, and P, is the usual Legendre function. The latter 
function can also be expressed in the usual way as a sum of terms involving 
the Legendre functions of cos 4, cos @ and cos v(¢, — @»2). For the 
ground state of helium and similar systems in |S states this is not necessary, 
as the distribution in three-dimensional space is spherically symmetrical 
and the wave function can only involve the radial coordinates and 6)». 
It is accordingly convenient to express the wave function in a form similar 


to (2.3), viz. 


v1, 2) = SSy,(1, 2)P,(cos 6,2). (2.4) 


When a function of this type is substituted in the wave equation (2.1) 
and use is made of the fact that a product of two Legendre functions can be 
expressed as a linear sum of Legendre functions, the wave equation sepa- 
rates out into terms each of which contains a Legendre function of the same 
order, say P,. Multiplying the equation by P, and integrating over the 
variable 4,2, all terms vanish except one set. In this way a set of coupled 
equations for the functions y,(1, 2) are obtained as follows 


(4 (#2) 452 (02)]+ 
oe » r, rs 
' 2Ln yw "On, rs Ore * One 


dF \ 


l 
+ | + Vin) + V(r) ¢ ¥, + 


ry ry 
Ect? Uy = Ey, (25) 


s, 8 


» 


where c\”’, are numerical coefficients, whose value is zero if r is negative or 
nonintegral. In particular the equation for Yo, the first term of the series 


(2.4) proves to be 


UT 707 240 LAO fig 0 
AF or (1 =) T r ar, (x 2)] + V(n) + Vin) Wo + 


 4U,/2r + Df ¥, = bY (2.6) 


r=0 


The Hartree solution for helium, or any method of representation of the 
(1s)* type, involves two approximations to (2.6). Firstly all terms in the 
expansion of 1/rj,. other than U» are neglected and secondly Yo (1,2) is as- 
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sumed to be expressible as a product of a function of the coordinates of elec 


tron | and the same function of the coordinates of electron 2. 

The calculations of Taylor and Parr* concerning the energy of the ground 
state of the helium atom, reported to this conference, are seen to fall within 
the above scheme. They express the wave function as a sum of terms, the 
first of which can be denoted symbolically by (1s, ls’), the second by a sum 
over 2p functions, the third by a sum over 3d functions and the fourth by a 
sum over 4f functions. These correspond to the first few terms of (2.4) ex- 
cept that the radical part of each term is expressed as a product of hydrogen- 
like functions with appropriate values of Z, the effective nuclear charge. 
In spite of this limitation the convergence of the calculated energy toward 
the observed one, as the wave function is represented by one, two, three and 
four terms, is remarkably good. The best value of the calculated total 
energy differs by only about one-eighth of one per cent from the observed 
value. Even a two-term expansion suffices to give an accuracy of about 
one-sixth of one per cent. 

From the wave function obtained by Taylor and Parr some information 
can be derived about the extent of spatial correlation in helium. This 
may be of two types. Firstly one electron may tend to be “‘out’’ when the 
in,’ and secondly one electron may tend to be on the opposite side 


‘ 


other is ‘ 
of the nucleus from the other. It appears from the first term of their wave 
function (of the 1s, ls’ type) that when 7, is small, the most probable 
value of r. is about 0.34 A, whereas when 1 is large, the most probable 
value of re is reduced to 0.23 A. This indicates the extent of in-out cor- 
relation for a function of this type. 

Some idea of the extent of angular correlation may be obtained from the 
first two terms of the wave function. When there is no angular correlation 
at all, the average value of the angle subtended by the electrons at the 
nucleus is 90°. The deviation from this value in helium does not seem to 
be large. It has a maximum when 7, and ry are both about 0.6 A and is 
then about 6°. When 7 and 7, have their most probable values (about 
(0.3 A) the deviation is about 4°. The average deviation for all values of 
r, and 72 is slightly less than this.‘ 

3. The Hydrogen Molecule.For the hydrogen molecule and similar 
homonuclear systems it is convenient to use elliptic coordinates. If &, and 
£ are the coordinates for the confocal ellipsoids, m; and 7 those for the hyper- 
boloids and $y the angle between two azimuthal planes, then 1/ry. can be 
expressed in terms of them by the Neumann expansion. This may be put 
in the form 


D 


c/2rn = OY D,,, (Ey &) Pim) Pim) cos vor (3.1) 
) 


r=Ov=( 


where c is the distance between the two centers and D,, is a function in- 
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volving only & and &. Each term is symmetrical for exchanges of the 
electrons and the symmetry properties of the various terms in the multidi- 
mensional space are determined by the factors in m, 2 and ¢@y. These may 
be symmetrical or antisvmmetrical for changes of sign of m or of m2, and 


there is periodic change of sign in dy. 
The wave function for the ground state of the hydrogen molecule, being 
a 'Z, state, must be symmetrical in the space coordinates of the electrons 


and so it is convenient to express it as a series of the same type as (3.1), viz. 


(1,2) = O (., + va,,) COS vor: (3.2) 
v=Q@ 
where y, , and y, , are each functions of the coordinates of electrons 1 and 
2, and have the same symmetry properties as the terms in (3.1), being 
either even in both m and m or odd in both n; and np». 

When y(1, 2) is substituted in the wave equation, the terms arising from 
(1/rijW(1, 2) can be arranged according to their symmetry types in the 
multidimensional space. All the terms which contain cos v@». can be 
equated separately to zero. Further the resulting set can again be sub- 
divided into those of symmetry y, , and those of symmetry y, ,. Since 
the expression must be identically zero for all points in the (£1, £, m, m2, d12) 
space and cannot change when m is changed to —m, the parts with y, , 
and y, , Symmetry must separately vanish. Thus a set of coupled equa- 
tions can be obtained for the set of functions y, , and y, ,, v ranging 
through all integral values. These have been given elsewhere by Len- 
nard-Jones and Pople.® 

Although direct solutions of these equations may not be possible at pres- 
ent, they form a useful background against which to examine the limitations 
of approximate functions. Thus the simple molecular orbital representa- 
tion of the wave function in the form f, 9(1)f, 9(2) involves the neglect of 
all parts of Y(1, 2) which depend on ¢,» and also of the part which has sym- 
metry Yo. It also expresses the wave function as a product and so auto- 
matically cuts out spatial correlation. A better approximation to the wave 
function would be given by 


y( i. 2) = fro ( | i 0 (2) My tuo ( | \fao ‘Gag (3.3) 


where /, 9 is a one-electron function independent of @ and antisymmetric in 
the plane » = 0. The negative sign in (3.3) is given because this proves 
to be necessary if f,, 9 is to be real everywhere. The condition that the 
energy should be a minimum leads to differential equations which f, » and 
fa o must satisfy.* 

Even the simple function (3.3) involves spatial correlation to some ex- 
tent. For if the position of electron | is fixed, the distribution of electrons 
depends on the square of a linear expression in f, 9(2) and f, o(2).. This 
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is not symmetric in the plane of symmetry and so there is a greater proba- 
bility of finding the electron nearer one nucleus than the other; in fact, 
the electrons tend to be on opposite sides of the plane of symmetry. By 
writing 


pF tla o) + ha) Vv y- xa = (fs - la) Vv 2 (3.4) 
the approximate form (3.3) is seen to be equivalent to 
W(1, 2) = xall)xn(2) + xe) xa(2). (3.5) 


The Heitler and London wave function and all similar functions involving 
directed atomic orbitals are seen to be particular cases of (3.5). 

A better approximation to (1, 2) can be obtained by ineluding in (3.3) 
terms 1n COS vd. ‘This leads to a generalization of (3.5) in the form 


¥(1, 2) = > 1XA> AL)xe, (2) + xa, .(Lxa, (2) feos vd. (3.6) 
, 0 
An expression of this kind implies not only correlation with respect to the 
plane of symmetry but correlation also in the azimuthal angle. It is still 
lacking in one respect. It does not allow for correlation of the “in-out” 
type, discussed above for helium. This would be included in the correct 
solution for ¥(1, 2) even in the first term of the expansion (3.2). 

Some indication of the effect of correlation on the energy of the hydrogen 
molecule for a wave function of the type (3.6) can be found from an analysis 
of some results recently published by Hirschfelder and Linnett (1950).? 
They use an approximate wave function built up from atomic orbitals of 
the hydrogen Is type, ¢4 and x, with adjustable screening constant. The 
first term, independent of the angle dy, is then 


dal(l)da(2) + dal(2)da(1) + ¥{ da(1)oa(2) + on(1)n(2)} (3.7) 


where y 1s an adjustable parameter. The second term of (3.6) is approxi 


mated to by taking x4; and yx, ; to be multiples of 


gira SIN 4; dary Sin Og 


which simulate atomic p functions. This part contains a numerical factor 
a which is adjusted to give the best value of the energy. 

By giving a and y various values, approximations can be obtained to the 
types of wave function dealt with above. Thus we have the following spe- 
cial cases: 

(i) Function 1 (a = 0,7 = 1): This corresponds to a product of simple 
molecular orbitals. 

(i) Function 2? (a 0, y adjustable): This is usually described as the 
electron pair function with tome terms. It is also of the form (3.3) pro 
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vided f, » and f, o are expressed in terms of atomic orbitals. The optimum 


value of y is 0.267. 

(iii) Function 3 (aand y adjustable): This takes in azimuthal terms in 
addition to those of Function 2. The optimum values of @ and y are 
—0.054 and 0.273, respectively. 

Using the tables given by Hirschfelder and Linnett,’ the various contri- 
butions to the total energy for these three functions can be calculated. 
These have been calculated for an effective nuclear charge 1.19 and a value 
of 1.47 for the internuclear distance, these being close to those given by 
them. ‘The results are shown in the table, the contributions to the energy 
of the molecule being expressed in kilocalories per mole relative to the 


separated atoms. 


CONTRIBUTIONS TO THE ENERGY OF THE HYDROGEN MOLECULE RELATIVE TO SEPARATED 
Aroms (IN KCAL. /MOoLEe) 
ELECTRON NUCLEAR #LECTRON- 
ELECTRON NUCLEAR NUCLEAR 
KINETIC POTENTIAL POTENTIAL POTENTIAL 
ENERGY BRNERGY ENERGY ENERGY TOTAI 
Function | 71 105 126 — 98] —79 
Function 2 87 379 126 — 985 —Q2 


Function 3 &7 372 126 —981 — 96 


The table shows up the improvement in the calculated binding energy as 
functions of different symmetries are introduced into the total wave func- 
tion. The second column shows how the average electrostatic energy of 
the electrons is progressively decreased, thus emphasizing the importance 
of spatial correlation. The first and last columns would be equal in mag 
nitude if the virial theorem were satisfied. The divergence is doubtless due 
to the fact that the values of the effective nuclear charge and the inter- 
atomic distance used were not exactly the optimum ones. 


* This paper was prepared in connection with a conference on ‘‘Quantum-Mechanieal 
Methods in Valence Theory” held at Shelter Island, September 7-10, 1951, under the 
auspices of the NATIONAL ACADEMY OF SCIENCES (see “Summary of the Conference”’ in 
a subsequent issue of these PROCEEDINGS) 

' Lennard-Jones, J. E., J. Chem. Phys. (in press) 

2 James, H. M., and Coolidge, A. S., [bid., 1, 825 (1933). 

§ Taylor, G. R., and Parr, R. G., these PROCEEDINGS, 38, 154-160 (1952) 

‘ A fuller discussion of the helium problem with particular reference to spatial correla 
tion is to be published elsewhere by the author and J. A. Pople (Phil. Mag., May, 1952) 

5 Lennard-Jones, J. E., and Pople, J. A., Proc. Roy. Soc., A, 210, 190 (1951) 

® The equations have been given elsewhere (Lennard-Jones and Pople, Proc. Roy. Soc., 
A, 210, 190 (1951)), where also a discussion of approximate methods of dealing with the 
hydrogen molecule is to be found 

? Hirschfelder, J. O., and Linnett, J. D., J. Chem. Phys., 18, 130 (1950) 
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GROWTH OF ANALYTIC FUNCTIONS ALONG A LINE 
By R. P. Boas, JRr.* 
DEPARTMENT OF MATHEMATICS, NORTHWESTERN UNIVERSITY 
Communicated by S. Bochner, March 29, 1952 


Let /(z) be regular and of exponential type less than in the half plane 
x 2 O, and let }A,} be an increasing positive sequence of unit density such 
that A,41 -— A, 26> 0. Itis known that under suitable restrictions and 
with a suitable measure of growth, the growth of the sequence {/(A,)} de 
termines the growth of /(x) for realx > 0. I have established the following 
results. 

Let e(x) be a concave increasing differentiable function such that e(x) = 
o(x), €'(x) S e(x)/x < Land log x = o(e(v)). Let tA,} satisfy the further 
condition that |A, — m S e(n). 

THEOREM |. /} 

7 log {(A,,) 
lim sup ; 
n> « €(A,,) 


then 


: log | f(x) 
lim sup 
png. ¥ ée(x) 


THEOREM 2. // 


: log! {(+A,) 
lim sup 
ome és €(A,,) 


where B is sufficiently large, then 
log | f(x) 


lim sup 


x 


e(x) 


THeoreM 3. If f(z) is entire and of zero exponential type, \f(+X,)} ts 
bounded, and €(x) satisfies the additional condition 


(x) Ss x ~2e(x) dx < 


then f(z) is a constant. 
THEOREM 4. /f f(z) ts entire and of zero exponential type and {( +X,) 
Olexp(AS)},0 <p < 1, then f(z) is of order not exceeding p. 


Theorem | interpolates between, but does not include, the theorem of V 
Bernstein! in which n An o(n) and the conclusion 1s 


° ° 1 . 
lim sup x! log f(x) lim sup A, log /f(\,)); 


n 
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and the theorem of Duffin and Schaeffer? in which ¢(x) is a constant, 1.e., 
< L, and the conclusion is that f(x) is bounded if [f(A,)} is 
bounded. Agmon® has treated similar rates of growth of the function but 
with a more restrictive condition on }A,}. As far as concerns the magnitude 
of e(x), Theorem 3 is an improvement on Levinson’s generalization‘ of 


n — X, 


Pélya’s theorem in which A, = n. Theorem 3 covers, for example, e(x) = 


x/tlog x (log log x)*} for large x, which is not covered by Levinson’s theo- 
rem, Levinson has shown that the convergence of (*) is essential. Levin- 
son's work is, however, more general in allowing complex 4, and a weaker 
spacing condition than |A, — A,,| 2 6,n # m. 

Theorems | and 2 have applications to the study of the behavior of power 
series on the circle of convergence, to Fourier series and to Tauberian 
theorems. 

The theorems are proved by a rather straightforward, but not entirely 
obvious, process® of comparing the given function with an interpolation 
series based on the points A,. This requires sufficiently precise estimates 
of an entire function with zeros at +A,, and similar functions, and it is these 
estimates which represent the essential novelty of the investigation re- 
ported here. The calculations are rather lengthy and will be published 


elsewhere. The essential estimates are as follows. Let 


« 


n l 


then (x + ty)| S Ae**™ for |y| < 1, and |@’(A,)| 2 Ae~4*%, where 
A denotes various positive constants. If e(x) satisfies the additional con- 
dition (*), then there is an entire function ¥(z) of arbitrarily small positive 
exponential type such that ¥(0) = | and log W(x) < -Ae(jx)). 

* Feliow of the John Simon Guggenheim Memorial Foundation. 

' Levinson, N., Gap and Density Theorems, New York, 1940, chap. VII. 

? Duffin, R. J., and Schaeffer, A. C., ‘Power Series with Bounded Coefficients,"’ Am. J. 
Math., 67, 141-154 (1945). 

§ Agmon, S., “Functions of Exponential Type in an Angle and Singularities of Taylor 
Series,”’ Trans. Am. Math. Soc., 70, 492-508 (1951) 

* Levinson, N., op. cit., chap. VIII. 

6 Cf. Boas, R. P., and Schaeffer, A. C., ‘A Theorem of Cartwright,”’ Duke Math. J., 
9, 879-883 (1942); Korevaar, J., ‘Approximation and Interpolation Applied to Entire 
Functions,"”’ Amsterdam thesis, 1949; Levin, B., ““O funktsiyah konechnoi stepeni, 
ogranichennyh na posledovatel’nosti tochek,’’ Doklady Akad. Nauk SSSR (N. S.), 65, 


265-268 (1949) 
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By LuTrHer P. EISENHART 
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Communicated March 17, 1952 


By definition a generalized Riemann space of coordinates x‘ (2 = 1, 
, n) is one with which is associated an asymmetric tensor g,,, whose de- 
terminant is not identically zero, and coefficients I“. of an affine connection 
asymmetric in the subscripts. It is understood that the coefficients I’, 
and I’,3 in coordinate systems x‘ and x’ are related by the equations 
O*x" , Ox or iy a 
= | (1) 


la '¢ + I Ve , ai af /. 
Ox “Ox ” "os “te” Oe 
In this equation, and throughout the paper the summation convention in- 
dicated by a repeated index, such as 7, j and y in (1) is used. 
Interchanging a and # in (1), we write the resulting equation in the form 
Ox" Ox’ Ox! i | Oenr 
yh = J? (2) 


, 


, , ji , 'g B a 
Ox %Ox'* Ox 7 Ox ' * Ox 7 


One-half the sum and one-half the difference of equations (1) and (2) 


are the respective equations 


a , Ox! Ox! oF 
J (3) 


P 8 aj ‘ — 4 ’ 
Ox’70x"! Ox’* Ox’! 2 Ox’? 


and 


, Ox! Ox’ » Oe 
y = 13 (4) 


ai ‘ _— ’ 
¥ Ox’* Ox’? ~ Ox’? 


/ . 
¥ 


yh a! . sph ly uJ af 
where [';; and I',3 are the symmetric parts of Ij, and [',3, and 'j, and [',} 


the skew-symmetric parts. This notation is used throughout this paper. 
Thus g;; and g,; are the symmetric and skew-symmetric parts of g,;.. From 


(4) it follows that I", are components of a tensor. 
The quantities 
Liss gn Vy - ginal jes (9) 
where g,,, , denotes the derivative of g,, with respect to x“ (a notation used 
throughout this paper) are the components of a tensor, being the covariant 
derivative of g,, in terms of the ['’s.'. Denoting by a,,, the components of a 
tensor to be specified later, we consider the equation 
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A A ' 
Sik * Lnjl th r Lul jk + ijk: ((}) 

From this equation we have 
vh wh ~ 
Liss s Lnil'sp + Lil, t yk (7) 


where d,, 1s the part of a,, symmetric in the first two indices, that is, 
| 


A jp = (Ajjp + jn). We choose a, so that 


vh vh 
C$ jx +- Lnjl 7 + Lal jk ~ 0, 


in which case equation (7) becomes 


sh sh 
Zirk = Sailing + Sia! 5 


from which it follows that? 
(10) 


hi . ‘ . - ° ° 
where |/,{ is the Christoffel symbol of the second kind in terms of the 
g,'s; and similarly for all the parts of the ['’s symmetric in the subscripts. 


In consequence of (8) equation (6) reduces to 
vh wh vii vh 
Py kt ™ Lrjl ik + Lil i. Lajl ik + g inl ik + tj 
where a,,, 1s the part of a,,, skew-symmetric in the first two indices. When 
equation (9) is subtracted from this equation, one obtains 
vi rh 
Lij “= Lng Vi + £ inl Ph + Cj jk (11) 

where g,;,, 1s the covariant derivative of g,; in terms of the Christoffel 
symbols, that is* 


h | | ‘ 
Ray /k Ligek T Bry tint + gin tie} (12) 


From equations (11) we have 


Valgisn + Bays + Bers) = BV in + alas + Oral Gi 
t+ /o(i yn H+ Ajng + eis). 
We replace this equation by the equations 
gayi + "/2 Gass + Cay ™ Bails + 1/2 Oya + Coy ™ 
1 


2 154i + “ii = 1/6(2 1; k + Zk /t + Lxi/j) (13) 


“anh 
Lux l ” + 


where the c's are components of a tensor skew-symmetric in the first two 
indices such 


Ciky + Ceyji + Chih ~ 0. (14) 
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Substituting from (14) in equation (11) we have in consequence of (14) 


Cae = . 3(28 45 / + L£ej i + Li — ijk r : 2( CL jx | + Axis). (15) 


These expressions satisfy (14). From (13) and (15) we obtain 


1 


vh , 
Ln Vi; Agia H+ Best + Susy) — Vol@ise + Get + Gexs). (16) 
If a, is chosen so that a,j, 0, we have 


yh 
Zn V Gj = | ALiyse HF Lagi tH Litsj)- 


If on the other hand we choose 


sh yh 
Gin = 2njVex + (Sin + BadT Sj, 


the condition (S) is satisfied, and 
fe vh vh 
QC ijk - Lil a Lil jk 
In this case (16) becomes 
vh 
LuV; =! 22 ij k + Lj i + Lu j)- 
If again we choose 
‘ sh 
ijn = 22 nl kis 
the condition (8) is satisfied, and in this case we have from (16) 
sh 1 


Lrxl a (Pi; k + Lj t + Lu joe 


When the condition of integrability 


oO ( O?x" )- O ( O?x" ) 
Ox’? \Ox’*0x"®J Ox’? \Ox’0x” 


of equations (1) is reduced by means of equations of the form (1), one finds 
ae wh . 
that the quantities I’,;, defined by 
sh = yh vh woah ado ph ‘ 
Dige = Tye — Ving + TET — Tiel y (21) 


are components of a tensor. In like manner from the condition of in- 


0 ( O%x" )- D) ( O%x" ) 
Ox’? \Ox’*0x"®F Ox \ox’70x"* 


of equation (1) one finds that the quantities TM, defined by 


tegrability 


alg _ yh yh vi ph wt ph 
Paik i ve is Vy.s + Tole - 9 


are components of a tensor. 
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Interchanging a and @ in the above conditions of integrability and apply- 
ing them to equation (2), one obtains the respective tensors 
hh wh hh wl qh 
"lise = ys a. Vik. i = A ee, ry Vie 
h al fica 
ijk = 1”, ji, k Ti j + I; Th =< Ppl j- 


In consequence of (10) equation (21) may be given the form 


vh. h vA sh yl pa 917 
be : - Riis + I a D ik /j > I if i = PyP ie (21 ) 


sh , . . . wah « ° a0. ° 
where I'/;,, is the covariant derivative of Ij; in terms of the Christoffel 


syinbols, and R¥,, is the Riemann tensor defined by* 
Rij. ss tae af tan + bie} Lit 
In like manner equations (22), (23), and (24) may be given the forms 
ie = —Ria +The —Tha +TGlh —TETh, 22’) 
Tye = —Riat Ve — lhe +P ule —THTE,  (23') 
? ‘ik on — Rix a Vii , oo Phi /3 7 Vin wi = Pil j- (24’) 
One-half the sum of the tensors (21) and (22) is the tensor 
Thin = Vie — /2Piey + Psd) + Pe — eT + PPh). 
Contracting for h and k, we have 
Ty = Tin — '/T i, + Ti) + Pe — PETG — Pig — Phe )- 
The expression in the second parentheses being a tensor, the other part of 
the right-hand member is a tensor; it is the Einstein tensor R,,.° 


From (21’) and (24’) one has the tensor 


ee rr" 


F (0 tik ca Ts.) = Dj k 


and from (22’) and (23’) the tensor 
: wh 
1/o( Tig nites " ijk) 


* See these PROCEEDINGS, 37, 311-315 (1951). 

' Eisenhart, L. P., ‘‘Non-Riemannian Geometry,” Am. Math. Soc. Colloquium Publica- 
tions, 5, p. 5 

*R.G.,p.17. This refers to the author's Riemannian Geometry, Princeton University 
Press. 

*R. G., p. 27. 

*R.G., p. 19. 

® Einstein, A., The Meaning of Relativity, 3rd ed. (revised), Princeton University 
Press, p. 114. 
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ON THE STATISTICAL THEORY OF TURBULENT DIFFUSION* 
By FRANGOIS N. FRENKIEL 


APPLIED Puysics LABORATORY, THE JOHNS HOPKINS UNIVERSITY, SILVER SPRING, 
MARYLAND 


Communicated by J. von Neumann, April 80, 1952 


1. Introduction.-The fundamental equations of turbulent diffusion 
have been first established by G. I. Taylor. In his paper on diffusion by 
continuous movements,' he considers the one-dimensional case of diffusion 
and finds the variance of the distance through which the particles have 
dispersed from a point-source. This variance is expressed as a function 
of the dispersion time and of two statistical characteristics of the turbulent 
field: the variance of the turbulent velocity and the Lagrangian correla- 
tion coefficient. Taylor's equations can be used as an approximation in 
the case of the three-dimensional turbulent diffusion from a continuous 
source in a fluid flow, when the turbulent velocities are small compared 
to the mean flow velocity.” However, in many applications to fluid 
dynamics, meteorology, oceanography, chemistry, ete., turbulent velocities 
are of the same order of magnitude as the mean velocity or even of a larger 
order. The main purpose of this communication is to consider the tur- 
bulent diffusion in a field of homogeneous and isotropic turbulence when 
the turbulent velocities cannot be assumed small as compared to the mean 
velocity. 

2. Fundamental [:quations of Turbulent Diffusion.—Let us consider 
within an incompressible turbulent fluid at rest, of homogeneous and 
isotropic turbulence, a fluid particle which at an initial time O is placed 
at the origin of the coordinate axes. Let u(t), v(t) and w(t) be the x-, 
y- and z-components of the instantaneous velocity of the particle. After 
a time /, the coordinates of the particle are 


; . . 
x= S u(a)da; y = fi'v(alda; 2 = f,'wlajda. 


The coordinates x, y and z are random variables’ and we shall investigate 
their variances as functions of the dispersion time ¢. The mean values of 
these coordinates Z = ¥ = 2 = O and their variances are functions of the 
dispersion time, and since the turbulence is homogeneous and isotropic, 
v= y? = 22 = Yl). 


We shall determine Y(t) as a function of the turbulence characteristics, 


We find 


y? = [fy 'v(a)da|? = So So 'v(a)v(a)day dar 


and introducing the Lagrangian correlation coeflicient' R,(h) 
v(tjo(t +h) /v*, 
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¥ = vf So'Ri(as oon a)day das. 
Taking into account that R,(h) is an even function 


y? = y? Sida fin" Rilada = v{fi'da fu“R(ada + 
Soda a "Ri (a)da). 
and after inverting the order of integration and determining the new limits 
of integration, we find 
y? = 2y? Sit — a)R,(a)da. (1) 


This relation was found by a different method by J. Kampé de Fériet® 
and is simpler for applications than the equation 


y? = 2y? Sc day Jo™ Ry(ar)day (2) 


originally demonstrated by G. I. Taylor.! 

We shall consider two special cases: when the interval of time ¢ is large 
and when it is small, as compared with the Lagrangian scale of turbulence 
Lp = Ss R,(a)da. 

Large Dispersion Time: When t > L, equation (1) becomes 


(t > Ly) y? = Qy? Lat — 20? fy” aR, (ada, (3) 


where the last term is constant. For very large dispersion times f, this 
constant becomes small as compared with the first term on the right 


and! 


(t> Ly) y? = 2y? Ly. (4) 


Small Dispersion Time: Consider now the case when t < Ly. Ex- 
panding the correlation function in a power series, taking into account its 
evenness, and neglecting terms of 4 of an order higher than the second, 
we find R,(h) ~ 1 — [h?/(2d,")] where d, is called the Lagrangian micro- 
seale of turbulence. Substituting R,(h) in (1), we obtain 


l #? re Z 
(i 1,) y= (1 — 12 -) vt*, (5) 


When the ratio //A, is small, the second term in the brackets may be 


_ neglected and! 
(t < Ly) y' 242, (6) 


General Case: When the dispersion time ¢ cannot be considered small or 
large as compared with the Lagrangian scale of turbulence, the variation 
of y* as a function of ¢ depends on the shape of the correlation curve. When 
this shape is known it may be possible to represent the correlation curve 
by a correlation function, Various representative functions® can be used 
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for this purpose and the variation of y* as a function of ¢ can then be found 
by using (1). 

3. Mean Concentration Distribution.—Let us consider a turbulent fluid 
at rest and let us assume that a large number of fluid particles are con 
centrated at an initial time zero at the origin of a set of rectangular axes. 
After an interval of time ¢ these particles are dispersed in a cloud by the 
influence of the turbulent fluctuations. If the turbulence is homogeneous 
and isotropic, the cloud of particles is spherically symmetrical. 

In the preceding paragraph we have given the variance y? for such a 
cloud of particles as a function of the characteristics of turbulence. We 
shall now investigate the distribution of particles within a cloud. This 
distribution can be expressed by Prob[x < a <x + dx, y < B<y + dy, 
s<y<s+dz| = ¥(x,y,2)dxdydz. The function P(x, y, 2) isa probability 


4a 
density function and f ff Y(x, y, z)dx dy dz = 1. We assume’ a 


three-dimensional Gaussian distribution and, therefore, a probability 
density function of the form 


3 


v(x, ¥, 2) = —, exp [—k*(x? + y? + 2%) ]. (7) 
v 


When the variance y? for this distribution is determined, we find k? = 2y?. 

4. Instantaneous Point-Source of Dispersion in a Turbulent Fluid at 
Rest._-Let us consider the case of a turbulent fluid at rest. At an initial 
moment 0, the particles are Concentrated at the origin of the axes. The 
number (or quantity) of partitles emitted at the initial moment is denoted 
by Qo. So(x, y, 2, t) is the mein concentration of particles at a point x, y 
z, and after a dispersion time ¢. The mean concentration is measured 
by the number (or quantity) of particles per unit of volume. We assume 
that Qy is very large (or that the particles are divisible) and, therefore, at 
the initial time §)(0, 0, 0,0) = ©, So(x, y, 2, 0) = 0. 

Applying (7), we find that after a dispersion time ¢, the number of 
particles found in a volume element dxdydz is equal to 5o(x, y, 2, t)dxdydz 
with 

; Qo eee ee ae 
So(x, ¥, 2, t) = (Oxy) exp | - er (x? + y? 4 ».| (S) 
In this equation y? is to be expressed as a function of v?, R,(h) and ¢ by (1). 

When we consider the dispersion of particles in a turbulent flow whose 
mean velocity 4 = UL’ is parallel to the x-axis, then the coordinate x in (8) 
should be replaced by (x — U1). 

5. Continuous Point-Seurce of Dispersion in a Fluid Flow.—Consider 
now, in a fluid flow, the dispersion from a continuous point-source emitting 
steadily Q particles per unit of time, The origin of a set of rectangular 
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axes 1s placed at the point-source with the x-axis parallel to the direction 
of the mean velocity 4 = LU. Due to the choice of axes, the components 
of the turbulent velocity are u’ = u — U,v’ = v, w’ = w. The mean 
concentration distribution for the time interval (¢, 4 + dto) of those Qdty 
particles which were emitted during the interval of time (t, f + dt) can 
be represented by (8), if we replace S)(x, y, 2, t), Qo and x, respectively, 
by S(x, y, 2, t — to)dto, Qdto and (x — Ut). Taking the sum of the effects 
of the emission from — © till /, we find that 


Su(x, y, 2) = + ae Six, y,2,t — a)da = S ” Six, y, 2, B)dBp (9) 


is independent of the dispersion time. §,(x, y, 2)U dydz is the average 
number of particles which pass per unit of time through an element of area 
dydz of a plane perpendicular to the direction of the mean velocity U. 
From (8) and (9), we find the mean concentration flux 
: ; ' Q ie I 5 — * - 
Su(x, y,2)U = l 7, exp { — - [((x — UB)? + y? + 27]? dB 
0 (2ry?) | 2y? f 
(10) 


where y? = 2y? So" (B — y)Rily)dy. 

Non-dimensional Mean Concentration Flux: We shall now represent 
(10) in a non-dimensional form. With this representation it will be easier 
to examine the dependence of the mean concentration flux on the various 
parameters. 

We introduce new coordinates & = x/UL,, 9 = y/ULy, & = 2/ULy. 
The Lagrangian correlation function is represented as a function of h/L, 
by &,(h/L,) = R,(h), and the intensity of turbulence is given by 7 = 
Vul?/U = Vy9"2/U = Vw"/U. 

The mean concentration flux is represented as a function of the above 
coordinates by a,,(&, n, ¢) = U°L,*[5y(x, ¥, 2)U/Q], and fas Grd a,,(&, 
n, ()dn dé = 1. With the non-dimensional notations, (10) becomes 


ro | l 
(t,9,¢) = eee - eps, - = x 
Sere (29) me 3( By) ( 27%2(B;) 


ie= ++ ep dai (1) 


where «?(8;) = 2 Jo*'(Bi — vWS.(ydy represents «?(r) for r = t/L, = By. 
We see immediately that the mean concentration distribution is a func- 
tion of (in addition to the coordinates): (a) the shape of the Lagrangian 
correlation curve determined by the function ®,(h/L,); (6) the intensity 
of turbulence 7 and (c) the product l’L, of the mean velocity by the 
Lagrangian scale of turbulence. 
We shall consider two cases: when the distance from the emission source 
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x is large and when it is small compared to L’L,, i.e., when ¢ is, respectively, 
large and small compared to unity. 

Large Distances from the Source: When the dispersion time is large, 
(4) can be applied and «? = 27. Replacing ¢°(3,) in (11) by 26), we find 


E r 

o,(& 7, ¢) = exp 
uw \s> ’ ore rs) 
n1™r 25: 


with r = Ve + 92 + ¢2. 
The dispersion factor « was determined neglecting the diffusion in the 


after integration 


direction of the mean flow. The value of « corresponds, therefore, to a 
one-dimensional case of diffusion. We now recompute c taking into ac- 
count the three-dimensionality of the dispersion wake. The new dis- 


persion factor t represents the standard-deviation \ ya” of the mean 


concentration measured in a plane perpendicular to the direction of the 
mean velocity along a diameter passing through. the dispersion axis Ox. 
We have 


1 ya? S * na (t, n, Ody 
L,? v? T°? Sf, “a, (&, n, O)dn 


Applying (12) and replacing in the two integrals the ratio »/£ by a hyper- 
bolic sine and using the Schlafli transformation, we find’ 

ta? = 2 es ae (13) 

Ko(&/2T?) 
Since the ratio of the first order Bessel function A, to the zero order Bessel 
function Ay tends to one when their common argument tends to infinity, 
we have 
lim tg? = 2é. (14) 

£/2T?—» a 

At the limit, we have, therefore, a relation corresponding to the equa 
tion which G. I. Taylor has given® for the one-dimensional case. 

It should be noted that equation (14) is valid only when & is much larger 
than unity. On the other hand, when £/27” is very large, equation (13) 
is well represented by its asymptotic relation (14). Equation (13) will, 
therefore, be most useful when at the same time both £ and 7 are large. 
More particularly, equation (13) should be used when the turbulent 
velocities are of the same (or of a larger) order of magnitude as the mean 
velocity. 

An approximate relation can be given for o,,(é, n, ¢) for the case when 
n and ¢ are small compared with & Developing (& — r) in a series and 
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neglecting the high order terms, we find’ from (12) 
tq 1 7? + r) 
¢Aat) = -——~ epi - — ~o 15 
im) © rae p ( aT? gt, — 


Small Distances from the Source: When the distance from the source 1s 
sufficiently small, we find from (6) that « = +r. After replacing in (11) 
i(f;) by 6, and pertorming the integration, we find? 


t - I | I lx g 
oulss m5) (27) +g “*P ~ oT? ie \2 Tr * 


1 DW ss 
exp O72 erfe {| — Vern (16) 


where r° £2 + 9? + ¢? and erfe (X) = 1 — (2/n) -* exp (—a?’)da. 

When the intensity of turbitence 7 is small compared to ({/r), the second 
term in the brackets is large com:pured with the first term and neglecting 
the first, we find 


Poet es g* | ¢ se 
eoXhs 5) * ae ren? | - 27? ( xe =) mm (- V2T r) ae 


For small values of » and ¢ as compared with € (and when 7 is small), 
we find the simple relation 


; i 1 vty? 
o,,(&, 1 $) exp ( — (18) 


QnT? ¢ 


oT? 


It should be noticed that only with equation (16) do we find finite mean 
concentrations upstream of the dispersion origin. 


* This work was done under Contract NOrd 7386 with the U.S. Navy Bureau of 
Ordnance. 

! Taylor, G. 1., Proc. Math. Soc. (London), 20, 196 (1921). 

Taylor, G. 1., Proc. Royal Soc. (London), 1§1-A, 421 (1935); Dryden, H. L., Ind 
Eng. Chem., 31, 416 (1939). 

3 A random function depends on a parameter w chosen at random in a measure space 
with a measure equal to 1. ‘The mean value of a function ¥(x, y, 2, ¢, w) is then given 
by ¥(x, y, 2, ¢) = SN (x, y, c, t, w)du. In physical applications it is not possible to 
measure such mean values and time or space averages are used instead. 

* The only correlation coefficient used here is the Lagrangian time correlation coeffi 
cient, found when the turbulent velocity of the same fluid particle is considered as a 
function of time 

® Kampé de Fériet, J., Ann. Soc, Sct. (Bruxelles), 59-1, 145 (1939). 

* Various representative correlation functions of the form exp|—&( {A} )"|¢(h) with 
n = | and 2 are discussed in a paper prepared by the author in 1942 and published as 
O.N.E.R.A. (Paris), Rep. Tech. 34 (1948). 

7 Frenkiel, F. N., U.S. Naval Ordnance Lab., NOLR., 1136, 67 (1949). 

* Roberts, O. F. T., has given a similar relation [cf. p. 36 of Sutton, O. G., Atmvos- 
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pheric Turbulence, Methuen (London) (1949)]|. To find his relation Roberts uses the 
differential equation of turbulent diffusion (OS; /Ot) + U(O0S)/Ox) = AV%py. It can 
be shown that the coefficient of eddy diffusion K is equal to v’?/,, if the molecular dif 
fusion is neglected. The differential equation is valid only when ¢ >> Ly. It is, how 
ever, possible [Frenkiel, F. N., Proc. Micrometeorological Symposium, International Union 
Geodesy and Geophysics, Brussels, Belgium, August, 1951] to define an “apparent co 
efficient of eddy diffusion” K, = (v’? jf — a)Rp( a) da, depending on the dispersion 


time ¢, which, when used instead of A in the above differential equation, makes the 


equation valid for all values of ¢ 


THE NUMBER OF SOLUTIONS OF CERTAIN EQUATIONS IN A 
FINITE FIELD 
By L. CaRLitz 
DEPARTMENT OF MATHEMATICS, DUKE UNIVERSITY 
Communicated by E. T. Bell, April 18, 1952 

Exact formulas for the number of solutions of certain special types of 
equations in a finite field have been obtained in recent papers by Hua and 
Vandiver,! Faircloth,? Faircloth and Vandiver;* see also Ward.‘ In the 
present note we state some additional theorems of this kind. 

1. Leta, ..., a, be rational integers such that (a), ..., a.) = 1; then 
we can construct a matrix of integers (a,,)), 1, 7 = 1, ..., Rk, such that 
ay = a;and |ay| = 1. Nowif Fis an arbitrary field, consider the equa- 
tion 


&%...&°!§ = a@ (a, & € F) (1.1) 


for fixed a. For a = 0 the solution is evident; for a # 0, no & = 0. 
We put 


m 
which is equivalent to 
ny ((bi3) = (ayy) 7). 


Thus (1.1) is equivalent to m, = a, 9, arbitrary but #0 for 7 = 2, ..., k. 


2. Now let F be the finite field GF(q), ¢ = p", and consider the equation 
net = f(b -. 2 Sr) (Eu Sy) € GF(Q)), (2.1) 


where f denotes an arbitrary polynomial with coefficients in GF(q). Ap 
plying the result of the previous paragraph we get 
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THEOREM |. Let N denote the number of solutions of (2.1) in &,, ¢,; let 
A(a) denote the number of solutions of 


os »r) =a (a € GF(q)); 
then 


N= (q' — (g = 1)*)A(O) + (g ~ 1)" 'S* A(a). (2.3) 


ari 


When an explicit formula for the number of solutions of (2.2) is available 
then (2.3) will furnish an explicit formula for the number of solutions of 
(2.1). In particular using the well-known results concerning 


OG; tf) = = (2.4) 
where Q denotes a quadratic form with coefficients in GF(g) and with 
discriminant 6 # 0, we obtain 

THEOREM 2. The number of solutions of 
&,"" &,"4 _ Ot, ’ Cy) (a; 2 1) 
forr = 2515 
er a A(q*(q' -g"') —-(¢- 1)*q*) 
where X = +1 or —1 according as (—1)*%6 is a square or a non-square of the 
field; forr = 2s + 1 the number of solutions ts q**** or gt *** + Ag'(q — 1)* 


according as (dy, , a,) ts odd or even. 
More generally we can determine the number of solutions of 


fe . te* wo OE, ...., 0) (a; > 1), (2.5) 


$l Sk 


in £,, U’, where the lL’, are polynomials e GF \q, x| of degree < m; for brevity 
we consider only the case r even. Using a result of Eckford Cohen® we 
obtain 

THEOREM 3. The number of solutions of (2.5), with r = 2s, is 


m 


1 
(q* =. (q ne 1") yams + (q° ma Ng” 1)(2m—1) a. 7g 2(8 2) +} 
z@=0 


(q — 1)*(q* — A)g@-POm-H, 


where \ has the same meaning as in Theorem 2. 
Similar results can be obtained if one or more of the LU’, in (2.5) are 
restricted to be primary of degree m. 
3. We consider next the equation 
‘, 
7 a, It," =a (a, a, £5, € GF(q)), 


rel j=l 


where 
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(da, . » Liz) ad | 


We have 
THEOREM 4. The number of solutions &,; of (3.1) ts 


gh! + ig IL (g® — gig — 1)**>4), (3.3) 
t= 1 


wherek =k +... tk, andl =q-—1 fora =0,t = —I1 fora # 0. 
In particular the number of solutions of 


a, a, b . 
angi 'm + + a,§;"'m' = a ((a,, 0:) = 1) (3.4) 


as gu) + ig. 
Let P(é) denote the left member of (3.1) and consider the equation 


P(&) as Hoy, ’ er (3.5) 
with f as in (2.1). Then we have the following extension of Theorem 1. 
THEOREM 5. The number of solutions of (3.5) is 


R(O)A(O) + RO) & Ala), 


ax 


where R(a) denotes the number of solutions of (3.1) and A(a) the number of 


solutions of (2.2). 
The second part of Theorem 4 leads to 
THEOREM 6. The number of solutions of 


P(é,n) = Qh, ...5 $7), (3.6) 
where P(t, n) denotes the left member of (3.4) and Q denotes a quadratic form 
is given by 


§q”** 1 4+ A(q — L)q* ** 1 


25-42 


\q? 
where \ has the same meaning as in Theorem 2. 

We can also obtain a result like Theorem 3 if in (3.6) we replace ¢, by 
U’;, polynomials in GF{q, x| of degree < m. However, we shall omit this 
generalization and state instead the following 

THEOREM 7. Let f denote an arbitrary polynomial with coefficients ¢€ GF(q) 
and let Q denote a quadratic form; then the number of solutions of 


= Of, ..., £2) (3.7) 


A(q?~' + X(q* — g’")) + (¢ — A)(q?-! — Aq’), 


where A denotes the number of solutions of f(t, ..., &) = O and X has its 
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usual meaning. In particular if the polynomial f never vanishes, then the 
number of solutions of (3.7) 1s 
ee ag) 
The corresponding result when Q is a form in 2s + 1 unknowns is more 
complicated. In the following special case, however, a simple result is 


obtained. 
THEOREM 8. T[f f 1s a polynomial that never vanishes, then the number of 


solutions , ¢, of 


P(é, ’ §;) ~ Ob, a | Cast 1) 


q'(q"* + mq’), 
where 4 = +1 or —1 according as (—1)*6 is or is not a square, and 6 ts the 
discriminant of Q. 
1. In the next place we state a theorem in which the condition (3.2) is 
weakened somewhat. We shall require some notation. For 8 ¢GF(q) put 


e(B) = e'"/? (8) = B+ BP+... + B”', 


G(B) = > e( pe?) (B # 0). 
é 


Also let (8) = +1, —1, 0 according as £8 is a square, a non-square or 
zero in GF(q). It follows that 


G(B) = ¥(B)G), GA) = qy(—1). (4.1) 
We now consider the equation 
ayes” + ....t+ a,t,""n,"* = Q, 


where (a,, 5,) = 2fort = 1, ...,7. 
THEOREM 9. The number of solutions &), ny of (4.2) is 


1+ q-! ze e(—aB) Il (¢ + (q — I)Pla,B)G(1)). (4.3) 
t=] 


BO 


By means of (4.1), (4.3) can be simplified; however, we shall not take 
the space to write down the final formulas. We remark that in place of 
(4.2) we can consider an equation of the form (3.1) with (aa, . . Gu.) = | 
or 2, but the final formulas for the number of solutions are rather com 
plicated. 

5. Theorem II of reference 1 suggests the following additional result. 

THEOREM 10. Let (au, .. dun) = d,, (dj, dj) = 1 fort # 7. Then 
the number of solutions of 
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with IN, ; &, A O, is 
ge "(q — 1)((¢ — 1)! + (-1)') (k = k, 


The total number of solutions of (5.1) ts qt *'~', where k 


' These PROCEEDINGS, 35, 481-487 (1949) 
2 Thid., 37, 619-622 (1951) 

8 Tbid., 38, 52-57 (1952). 

4 Thid., 37, 113-114 (1951) 

5 Duke Math. J., 14, 543-557 (1947). 


ON SPECTRAL ANALYSIS OF STATIONARY TIME SERIES* 
By ULF GRENANDER AND MURRAY ROSENBLATT 


DEPARTMENT OF MATHEMATICAL STATISTICS, UNIVERSITY OF STOCKHOLM, AND CoM 
MITTEE ON STATISTICS, UNIVERSITY OF CHICAGO 


Communicated by Saunders MacLane, April 30, 1952 

1. The present statistical theory of analysis of stationary time series 
(e.g., extrapolation) has assumed complete knowledge of the covariance 
sequence or equivalently of the spectrum of the process. It is, therefore, 
of great importance to be able to estimate one of these. However, knowl- 
edge of the spectrum seems to yield greater immediate insight into the 
structure of the process. This seems to have first been noted in a funda- 
mental paper by Bartlett.'. An unpublished paper by Tukey® deals with 
some aspects of the problem of estimating the spectrum. 

2. Consider a discrete parameter stochastic process x, (f = ..., —1, 
0, 1, ...) of the following form: 


where }a,} is a real sequence and 


PS ae 


x 


The £,’s are assumed to be independent and identically distributed real- 


valued random variables with 
Et, = 0, Et,? = o?, EE = py. 


The conditions cited above imply that the stochastic process x, is strictly 
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stationary. Such a process is known to have an absolutely continuous 
spectral distribution function F(A) with a density given by 


| a. 
F’(X) =f(Xd) = - pF ae a’. 
a x 


Proofs of the results cited below will require f(A) to satisfy a smoothness 
condition and the &,'s to have moments up to the eighth order. 

As the process is real-valued, the spectral density is symmetric about 
zero and hence one need only consider a spectral distribution function 
F(X) as the integral of the spectral density from zero toA, O< A< am. A 
reasonable estimate of the spectral distribution function F(A) is 


Fy*(s) = So* In(Ddl 


where 


| 
Iv(l) = on 


is the so-called periodogram. Professor Doob suggested that we investi- 
gate the problem of giving confidence bands for the estimation of F(A) in 
terms of Fy*(A\). In the case of a normally distributed process, Theorem 
3 yields a solution which is non-parametric with respect to the spectrum. 
In the general case, Theorem | and an appropriate estimate of e = “ seal. 
0 

can be used to obtain conservative bands non-parametric with respect to 
the spectrum. Theorem 4 can be used to test whether two independently 
drawn time series come from processes having the same spectrum. 

3. The theorems are as follows: 

THEOREM lI. 


. | . . | |- | 
lim P' max.| N[Fy*(\) — F(A)]| < al = P35 max. |Z(A)/< ! 
N—e (O0S)XSn \ (O<\Sx, \ 
where Z(X) is a normal process whose distribution is completely determined by 
the mean value E:Z(X) = O and the covariance 


EZ(\)Z(u) = eF(A)F(u) + 2G (min. (A, u)), 


GA) = Sfo* f*(Ddl. 
Corotiary. /fe = 0 


lim PS max.|//N[Fy*(A) — F(A)]| < al = P'S max. |W(A) < al 
| \ \ 


N—>o (O<SASe (05 <28G(n)! 


where W(X) 1s the Wiener process, 1.e., 





VoL. 38, 1952 MATHEMATICS: GRENANDER AND ROSENBLATT 


EW(A) = 0 EW(\)W(u) = min. (A, w) 


and W(X) ts normally distributed. 
Note that e = O when the &, are normally distributed. 
THEOREM 2. The statistic 
l kN@ 9 
Ry = ”. 
2 = —kna N?2 
converges to 29G(m) in probability as N—~ © ifk >0,0<a< 1. The 
c,'s are the product sums 


v 


N v| 
— 
€, si z XnXn lye 
1 


THEOREM 3. Ife = 0 


lim P\ max. |} N[Fy*(A) — F(A) ] } P 
JOSA<s - _\ = P) max.|W(A)!< atl. 
\ 


Na i a a < « 
| Ry 2 \ (OS\AS1 

It is well known that the probability on the right can be explicitly 
evaluated in a convenient form. 

THEOREM 4. Jf e = Oand if Fyi*(d), Fyo*(A), Rui, Ryo are the statistics 
defined above corresponding to two independently drawn samples from the 
same stochastic process, then 

lim P( max. VN[Fui*(\) — Fre*(a)] } : : 

° JO<ASx < a> = PP) max. WA) |S 
: (Ry + Ryo)" =) eee 


4. Although the plausibility of the basic Theorem | is apparent by 
heuristic arguments, the proof requires a long and formidable sequence of 
lemmas. The practical value of the theory cited above can only be 
judged after having completed an extensive computational program in- 
volving analyses of artificially generated time series. Some preliminary 
analyses have already been made and tables of the related distributions 
are being compiled. Considerable investigation must be made of the 
practical application of the theory. The proofs of the above theorems and 
considerable elaboration of the theory (e.g., solution of the corresponding 
problem for the spectral density) as well as aspects of the practical applica- 
tion of the theory will be contained in a forthcoming publication. 


* This work was supported by a grant from the Office of Naval Research. 

! Bartlett, M. S., ‘“‘Periodogram Analysis and Continuous Spectra,’’ Biometrika, 37, 
(1950). 

2 Tukey, J. W., “Measuring Noise Color,’ unpublished manuscript. 
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ON THE THEOREM OF RIEMANN-ROCH FOR ADJOINT SYSTEMS 
ON KAHLERIAN VARIETIES 


By KUNIHIKO KODAIRA 
INSTITUTE FOR ADVANCED StTuDy, PRINCETON, N. J. 
Communicated by Hermann Weyl, March 10, 1952 


1. Let Yt be a compact Kahlerian variety of (complex) dimension m. 
Suppose a (possibly reducible) compact analytic subvariety S of YW of 
dimension m — | as given and denote by YW (.S; Yt) the linear space con- 
sisting of all meromorphic m-ple differentials W on Yt which are multiples 
of —S in the sense that (W) + S 2 0, (W) being the divisor of W. The 
problem concerning the theorem of Riemann-Roch for the adjoint system' 
of S consists in expressing the dimension of the linear space W(S; Yt) in 
terms of “numerical characteristics’ of Yt and S. The purpose of the 
present short note is to give an answer to this problem? in the following 
two simple cases: (a) S is an irreducible non-singular variety, (b) S consists 
of two irreducible components S,, S, which are both non-singular and the 
intersection V = S; Sy is an irreducible non-singular variety having the 
intersection-multiplicity 1 in the sense that at each point qeI’ the two 
varieties S,, S) have different tangential planes. 

2. We start with the following remark: Let V be a compact analytic 
non-singular subvariety of Yt of (complex) dimension m — h and t be an 
arbitrary pure p-current® of type s on V. Then we can associate with t 
a pure (p + 2h)-current 7(t) of type (s + h) by setting 7(t)[®] = t[&y] 
for a “variable” (2m — p — 2h)-form ® of Class C” on Wt, where 9%, 
denotes the (2m — p — 2h)-form on V induced by &. Clearly the corre- 
spondence t — 7(t) is linear and one-to-one. Moreover we have d7(t) = 
T(dt), since (db)y = d(%y). In view of these facts, we identify t with 
7(t) and consider t = T(t) as a pure (p + 2h)-current of type s + h 
on IN. 

3. First we consider the second case (b). Every m-ple differential 
W € WS; MW) can be considered as a pure m-current of type 0 on W, since 
the integral W[¥] = ,QW-W converges absolutely for an arbitrary 
m-form W of Class C” on Yt. Similarly, every (m — 1)-ple differential 
w, € WI; S,) is a pure (m — 1)-current of type 0 on S, (A = | or 2) and 
consequently w, can be considered as a pure (m + 1)-current of type 1 on W. 
The exterior derivative dW of the current W has the representation dW = 
w + we, where w, € WT; S,) (A = 1, 2). This can be proved as follows: 
First, let q be a point on = S; 9 Sy and (24), 292, . . . 5 2gm) be a system 
of local coordinates on YQ with the center q such that S, or S, coincides 
in a neighborhood U(q) of q with the coordinate plane z,, = 0 or 2. = 0, 
respectively. Then W has in U(q) the expression Wo = (24,29) 
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Fd2q\d2q2 dZm, Where F, is a holomorphic function defined in U(q). By 
a simple calculation, we get 


dW(®) = 21}, 52 2| Fyd2xdz; -. dn ® 


So a 
- S,2 1 Fydaid2zz . . . d2m*®s,} 


Si 


for an arbitrary (m — 1)-form # of Class C” whose carrier lies in U(q), 


where 2 = 2, (k = 1,2, ...,m). Consequently, setting 
s , 1 ‘y o ° o o ° ° 
Wig = 2r1 2 ae F,(0, Sy2> Saas - -» Sam Jd Zqxd2q3 ; BS, 
‘ - i= 
oq = 2wt tq, Fy(2q1, 0, 2qa, . . «> Sqm) OZqidZq3 . . . AZqm 


and considering the (m — 1)-ple differentials w,, on S, % U(q) (A = 1, 2) 
thus defined as pure (m + 1)-currents of Type | in Ul(q), we obtain 


dW = arg + we, in U(q) (qeT). (1) 


Second, take a point p eS, — I (A = 1 or 2) and choose a system of local 
coordinates (25), Zp2, .. ., 2pm) With the center p such that .S, coincides in 
a neighborhood U(p) of p with the coordinate plane z,, = 0. Then W has 
in U(p) the expression W = 2);' Fydzyid2y ... day, where F, is a holo- 
morphic function defined in U(p). Moreover, setting wy = 2ri F, 
(0, 2,2, Zym) *d2y2 dZym, we get, in the same manner as above, 


dW=a4,, inl(p) (pe S, —TP). (2) 


It follows from (1) and (2) that dW has the expression dW’ = w; + we, 
w, € WT; Sy). 

Now we shall show that, for given differentials w, « WT; Sy) (A = 1, 2), 
there exists a differential W ¢€ XW(S; W) with dW = w, + w. if and only if 
w1, we Satisfy dw, + dw. = 0 and Ha, + Hw. = 0 as (m + 1)-currents on 
Wt, where Iu, denotes the harmonic part* of the current w,. It is obvious 
that the existence of W with dW = a, + we implies d(w: + aw) = 0 and 
IT(a, + w.) = 0. Assume conversely that d(a@, + w2) = H(w, + w) = 0. 
Then with the help of Green’s operator G and the Hodge-Weil’s operator A, 
we can construct an m-current® W = (—idA + 6)G(a; + w2) of Type 0 
satisfying dW = a, + w. Moreover this current W is a holomorphic 
m-ple differential in Yt — S. In order to show that W is a meromorphic 
differential belonging to YW (S; We), let p, 2%, 4, 24 have the same meaning 
as above. Then w has on S$, 9 U(p) the expression a, = /,(2p2, » Sm) 
d2Zy. ... d2ym, where f, is a holomorphic function defined on S, % U(p). 
Setting 6, = (24i)~'2))'fy(ap, - . -5 Spm)dZyrd2yo . . . d2ym, we get therefore 
a pure m-current 6, of Type 0 defined in U(p) satisfying d0, = w,. The 
difference W — 4, satisfies d( W — @,) = 0 and consequently W — 6, must 
be a holomorphic differential’ in U(p). This proves that 2,,W 7s a holo- 
morphic differential in \W(p) for each pe S, — V. For q eT, it follows from 
the above result that 242,,W is holomorphic in U(q) — Tl; hence, by 
virtue of a theorem of Hartogs,’ 2:,22,W must be a holomorphic differential 
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in \l(q) for eachqeV. Thus we see that W is a meromorphic differential 
belonging to WS; We). 

1. We denote by (I) the linear space consisting of all p-ple differentials 
of the first kind on I and set g,(M) = dim, (MM). In the first case (a), 
it can be shown in the same manner as above that, for every W ¢ WS; Wd), 
w = dW is an element of G,,_,(S) and that, for given w e ,,—,(S), there 
exists We WS, I) with dW = w if and only if Hw = 0, where each 
element of @,,—,(.S) is to be regarded as a pure (m+ 1)-current of Type | 
on 9. Let {B,)v» = 1, 2, » £m— (V0 be a base of the linear space 
(%,, ,(M) which is normalized in the sense that Sv B,:*B, = 6,,. Then 
Hw has the expression Hw =(—1)"~' >> w[B,]*B,. Consequently, there 


exists W € WS, Mt) with dW = w tf and only if w satisfies 
w(B] = 0, for all B €®,,— (M0, (3) 


where w is a given element of ,,_,(S). Thus, denoting by $(.S) the sub- 
space of ,,,(S) consisting of all w e G,, —,(S) satisfying (3), we infer that 
the linear mapping W—~>w = dW maps W(S, Wt) onto H(S), while dW 
vanishes if and only if W belongs to the subspace @,,(QM) of WS, M). 
Hence we get the relation dim W(S, IW) = g»(Wt) + dim H(S). On the 
other hand, w[B] = [By] = (—1)"~2%>??)""! fw-*Bs, and 
Bs ¢ ®,,—,(S) for all B eG, (MY. Hence (S) is the subspace of G,, — ,(.S) 
which is orthogonal to the subspace {Bs|Be@,,_,(M)}. Setting / = 
dim {B| Be @,,_ (M9, Bs = O}, we infer therefore that dim 6(S) = 
Zm—10S) — gm—-\(QVO +1. Combining this with the above result, we 
obtain the following 
THEOREM |. Jn the first case (a), we have 


dim WS, WM) = gm(M) — gm— 1M) + gm—10S) + 1, (4) 


where | is the number of linearly independent (m — 1)-ple differentials B of 
the first kind on YN which vanish on S in the sense that Bs = 0, Bs being the 
differential induced on S by B. 

Remark: In case w ¢,,_,(.S) satisfies (3), the current W = (—idA + 
5)Gw satisfies not only dW = w but //W = 0 and therefore W[P] = 0 
for all P ¢ G,, (M0. 

5. We turn to the second case (b). Let w° be the subspace of YW (1; S;) 
+ WP; S.) consisting of all a, + w., w, « WT; S,) (A = 1, 2), satisfying 
the conditions dw, + dw, = 0, Hw, + Hw, = 0. Then the result of Sec- 
tion 3 shows that the linear mapping W — a; + w = dW maps WS; Wd) 
onto w°, while dW vanishes if and only if W ¢,,(9). Hence we get 


dim WS; WM) = g,,(M) + dim w°*. (5) 


Clearly the condition //a, + Hw. = 0 is equivalent to the following one: 
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wi [B] + w[B] = 0, for all B e G,, _,(M). (6) 
Now, applying the argument exponded in Section 4 to each pair of varieties 
S,, PT (A = 1, 2) instead of IY, S, we infer that, for an arbitrary element 
w; + we ew’, the exterior derivative & = dw, = —dw, belongs to G,,_2(T) 
and satisfies as a current on S, the relation 
Ela] = 0, for all ay eG, —9(.S,) = y (7) 
and that, if an element £ of G,,_(I") satisfies (7), there exist w, « W (Tl; S,) 
(A = 1, 2) satisfying dw, = —dw. = § and 
w,[B,] = 0, for all B, « G,, — (S,) (A = 1, 2) (8) 
(see the remark at the end of Section 4). Obviously (S) implies (6). 
Hence the linear mapping w; + w.— £ = dw, = —dw. maps Ww” onto the 
subspace (I) of G,,—2(I') consisting of all & €@,, (I) satisfying (7). 
The kernel f of this mapping is given by 
f = lw aa We w, €Y,, (Sy), ol Bs | + w2[Bs, | = 0 (Be G,, (WE) }. 
Setting / = dim |B) Be G,,— (M0), Bs = Bs, = 0}, we get dim f =g,,_, 
(S1) + £m—1(S2) — gm— (MN) + 1 and therefore 
dim w° = gy, —1(S1) + om—1(S2) — gm— (M0) +2+ dim H(T). (9) 
It is obvious that the condition (7) 1s equivalent to 
Elaup ia aor | aa 0, for all a € (S,, o(.S)) (A = & 2). (10) 
Hence, setting g»—o(S) = dim fay + ay! ay, €G,,—2(Sr), arp = aap}, we get 
dim D(T) = gm—o(T) — 2m—2(S1) — gm—2(S2) + gm—o(S). 


Combining this with (9) and inserting it in (5), we obtain the following 
THEOREM 2. In the second case (b), we have 


dim WW(S; M) = gua(M) — gm—1(DO) + gm—1(S1) — 2m —2(S1) 
+ teal) ~ te fd * es eee 


where | is the number of linearly independent differentials B eG, — \(M) 
satisfying Bs, = Bs, = 0 and gm_2(S) is the dimension of the linear space 
fay + ay! a €G,—2(S,) (A = 1, 2), arr = ae}. 

6. In this section we shall consider the case in which Yt is an (irreducible 


non-singular) algebraic variety imbedded in a projective space. Let / be 
an arbitrary hyperplane section of YW which is irreducible and non-singular. 
Then a theorem of Lefschetz’ asserts that, in case p S m — | (m = dim Wi), 
every p-cycle Z on Wis homologous to a p-cycle Zp lying on / and that, 
in case p S m — 2, a p-cycle Z» lying on FE is homologous to zero on E 
if Zz is homologous to zero on YW. It follows from this theorem that 
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(I) in case p S m — |, a p-ple differential P € (IN) vanishes identically if 
P,, vanishes identically on EF: and that (11) in case p S m — 2, the mapping 
P— Py, maps (I) isomorphically onto (12). The second proposition 
(II) implies obviously that g,(/1) = g,(IN) for p S m — 2. Incidentally 
we denote by K the canonical divisor on WN. 
Now, for an arbitrary irreducible non-singular algebraic variety V of 
dimension n, we set 
a(V) = g,(V) — gn-1(V) + Ba-20(V) — ... # (—1)"7*g,(V). (12) 
Then we have 
THEeoreM 3. The dimension of the complete linear system | K + E) on Wt 
is given by 
dim |K + E| = a(M) + a(F) — 1. (13) 
Proof: Since dim |K + E| = dim W(E; M) — 1, we get, by Theorem 
l, 
dim |K + FE] = gn (QM) — gm—\(M) + gm— CE) +1 — 1. (14) 


Now, it follows from the above proposition (1) that / = 0, while g,(/) = 
g,(W) for p Sm — 2. Hence we obtain from (14) the formula (13), 


q.e. d. 

THeoreM 4. Let S be an arbitrary irreducible non-singular subvariety 
of IW of dimension m — | and E be an irreducible non-singular hyperplane 
section of IR which cuts out on S an irreducible non-singular subvariety ES 
with the intersection multiplicity 1. Then the dimension of the complete 
linear system |K + S + E| is given by 

dim |K + §+ E| = a(t) + a(S) + a(k) + a(S) — 1. (15) 

Proof: By virtue of Theorem 2, we get 
dim |K +S+E| = Em(D) — gm— (DO + vm—10S) — g2m—2(S) 

+ Lm (2) = Lm a(t) + £m o( ES) + £m (S + E) + l aia if 

Now, we have g,(/) = g,(Q) for p S m — 2 and, since ES is a hyper- 

plane section of S, g,(/2S) = g,(S) for p S m — 3, while it follows from 
the proposition (I) that / = 0. Hence we obtain 

dim |K + S+ E| = a(M) + a(S) + a(E) + a(ES) — 1 + gn—o(S + E) 

— £n—2(D). 

It is sufficient therefore to show that) g,—o(S + FE) = gm—o(M). 
£m—2(S + £) is the dimension of the linear space 

W = fay + ay] are G,9(S), a2 € Gy AE), ay = any}, 


where [= /S. For an arbitrary element A ¢€@,,~.(M0), we have 
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Ase, (S), Ap eG, (FE) and (As)p = (Ag)p; thus A + Ay + Ag is 
a linear mapping of ,, (M2) into ®. It follows from the proposition (II) 
that, for every element a; + a, of (, there exists one and only one element 
A of G, (M0) with Ap = ae. This element A satisfies (a, — As)p = 
= ES is a hyperplane section of S. Hence, by 


’ 


doy — (Ag) = O, while I 
the proposition (I), we get As = a. Thus we see that A + As + Ax is 
a one-to-one mapping of @,,—.(Y) onto . Consequently we obtain 
Zm—2(M) = gm-2(S + E), g.e. d. 


In case there exists on M® at least one non-trivial meromorphic m-ple differential 
W», the canonical divisor K on YX may be defined by K = (Ws). Then the dimension 
of WS; IW) is equal to the dimension of the adjoint system |K + S| of S increased by 1. 

2 For the special case in which W is a Kahlerian surface, see K. Kodaira, ‘The 
Theorem of Riemann-Roch on Compact Analytic Surfaces,"’ Am. J. Math., 73, 813-874 
(1951), §6. 

3 See Kodaira, loc. cit., §1, where the reader will find a summary of the results in the 
theory of harmonic integrals which are necessary for the present note. 

‘de Rham, G., and Kodaira, K., “Harmonic Integrals’’ (mimeographed notes), 
Princeton (1950), pp. 63-68; see also Kodaira, loc. cit., p. 817. 

5 Kodaira, loc. cit., Theorem 1.4. 

6 Kodaira, loc. cit., Theorem 1.3. 

7 Hartogs, F., Acta Math., 32, 57-79 (1909) 

5 Lefschetz, S., L’Analysis situs et la géométrie algébrique, Paris, 1924, pp. 89-91. 

® It can be shown that a( V) is equal to the arithmetic genus of V. See K. Kodaira, 


nov or 


“Arithmetic Genera of Algebraic Varieties,’’ these PROCEEDINGS, 38, 527-533 (1952). 


ARITHMETIC GENERA OF ALGEBRAIC VARIETIES 
By KUNIHIKO KODAIRA 
INSTITUTE FOR ADVANCED STuDY, PRINCETON, N. J. 
Communicated by Hermann Weyl, March 10, 1952 
Let us consider an arbitrary irreducible non-singular algebraic variety 
V, of dimension n imbedded in a projective space, and denote by g,(V,) 
(1 S k Sn) the number of linearly ifidependent k-ple differentials of the 
first kind attached to V,. In the present short note, we shall prove a 
conjecture of F. Severi! to the effect that the arithmetic genus py(V,) of Vn 
is represented in the form 


Pal Vn) -_ £ nl V,) — £a 1( Vn) + £n 2( V ») ens 7 1)" "gy ( Fob 
(1) 


It follows from (1) that the arithmetic genus p,(V,) is a birational invariant 
of V,, since each g,(V,) is known to be a birational invariant? of V,. 
1. We set 
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r 


a(V,) “ gal V,) — £n (V,) T £0 ( Vy) 


Then we have to prove the identity 
pal V2) = a(V,). 


First we consider an arbitrary but fixed irreducible non-singular variety 
WM = V,, of dimension m = 2. A linear system | A} of effective divisors 
A on Yt without fixed components corresponds to a rational mapping?® 


~ 


of Wi into a projective space J, and this correspondence is one-to-one. 
We shall say that the linear system \A} is sufficiently ample if the corre- 
sponding rational mapping induces a biregular imbedding of YR into S. 
By virtue of theorems of Bertini, a general member A of a sufficiently ample 
linear system }A} is an irreducible non-singular subvariety of I. It is 
obvious that, if {A} is sufficiently ample, the complete linear system | A| is 
also sufficiently ample. Again, if |A} is sufficiently ample and if | B{ has 
neither fixed components nor base points, then the complete linear system 
A+B is also sufficiently ample. Assuming that J is imbedded in a 
(fixed) projective space S, of dimension d, we denote by F, the section 
of Yo cut out by a general hypersurface &, of order h. Then it follows 
from the latter proposition that, for a given divisor D, the complete linear 
system |D + E,| is sufficiently ample if h is sufficiently large. In what 
follows we denote by S, 7, UL irreducible non-singular subvarieties of Wt 
of dimension m — | and by AK the canonical divisor on YW. Incidentally, 
for an arbitrary divisor D on an irreducible non-singular variety V, we 
denote by Dy the complete linear system on V consisting of all effective 
divisors which are linearly equivalent to D. 

TuroreMm 1.4 /f S belongs to a sufficiently ample linear system on WN, 


then we have 
dim | AK + Sig, = a(Mt) + a(S) — 1. (4) 


THEOREM 2.° Assume that m = dim Yt 2 3. If S belongs to a suffi- 
ciently ample linear system on IN and if the intersection S-T is an irreducible 
non-singular variety with the intersection multiplicity 1, then we have 


dim | K + S+ 7} gy, = a(t) + a(S) + a(T) + a(S:T) — 1. (5) 


The above two theorems can be proved with the help of the theory of 
harmonic integrals.® 

THEOREM 3.) If | S| g; = | T| gn, then a(S) = a(T). 

Proof: In case m(= dim Wt) = 2, the genus a(S) of the curve S on the 
surface Yt is determined by the relation’ 2a(S) — 2 = 1(S,AK + 5S), 
where /(.S, A +S) denotes the intersection number of S and AK + S. 
Hence a(S) = a(7) if | Sig, = | T)q. Now, assume that our theorem is 
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true for m <n and consider the case in which m= n> 2. Let E = F, 
be a general hyperplane section of J. Then S or 7 cuts out on £ an irre- 
ducible non-singular subvariety S-- or 7-£, respectively. Since 
K+ S+E\y =|K+7+ Ely, we get, using (5), a(S) + a(S-E) = 
a(T’) + a(T-E). On the other hand, it follows from | S| g, = | 7'| g that 
|\S-Ely = |T-El|,. Hence, by the induction hypothesis, we get a(S: /) = 
a(7T’- FE). Consequently we obtain a(S) = a(7’), q. e. d. 

THeoreM 4. Let lL’ be an irreducible non-singular subvariety of MN of 
dimension m—1. If |S\yy is sufficiently ample, then |K +S + U| 
cuts out on LU’ a complete linear system. 

Proof: Denote by |K + S+ Uy: U the linear system cut out by 
(K + $+ U|\q on U and by Ky the canonical divisor on UL’. Clearly, a 
general member S of |S) gy cuts out on Ul’ an irreducible non-singular 
variety S-U with the intersection-multiplicity 1, and |S-U| y is suffi- 
ciently ample. Now, we infer readily that | A + S + U | m-WU is con- 
tained in the complete linear system Ay + S- lv, while 


dim {|K + S+ U|\m-U} = dim |K + 5+ U|g — dim |K + S|q- 1. 
By virtue of (4) and (5), we get 


dim (Kv +S§-Uly =a(U) +a(S-U) —1, 
dim | K + S| gp = a(t) + a(S) — 1, 
dim | K + S + Ug = a(M) + a(S) + a(U) + a(S-U) — 1. 


Consequently, we see that dim } [K+ 5+ Ulm: U} = dim (Ky + §-U| v. 
This shows that |AK + S+ Ul’ y-U coincides with’ the complete linear 
system |Ky + S-U) v, q. e. d. 

Coro.iary. Let LU be an irreducible non-singular subvariety of IN of 
dimension m — 1. Then, for an arbitrary divisor D on IN, |D + En\ gy 
cuts out on U a complete linear system if h is sufficiently large. 

2. We denote by /* the intersection of k general hyperplane sections 


of @. For each k < m — 1, F* is an irreducible non-singular subvariety 
of M of dimension m — k and, as one readily infers from Theorem 3, 
a(E*) is independent of the choice of £*, i.e., a(/:") is determined uniquely 
by IM, |) gq and k. EH” is a set consisting of a finite number of points 
on IR. We denote by a(k™) the number of points of E" decreased by 1. 
Clearly a(E”) + | is equal to the order of IM. 

LemMMA |. Assume that the intersection S-T 1s an irreducible non-singular 
variety with the intersection multiplicity 1. If | U\q = |S + T7|\m and if 
|S! qn and | LU! are both sufficiently ample, then we have 


a(U) = a(S) + a(T) + a(S-T). (6) 


Proof: The lemma is an immediate consequence of (4) and (5). 
THEOREM 5. We have 
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m 


a(E,) + (-1)"" = > (1) fa(E*) + (—1)™~-*]. (7) 
l 


h 


Proof: We prove the relation (7) by induction on m and h. In case 
m = 2, (7) is an immediate consequence of the relation 2a(,) — 2 = 
I(k,, ky, + K). Now, assume that (7) is true for m S n — | and con- 
sider the case in which m = n> 2. For h = 1, (7) is valid trivially. 
Assume therefore that (7) is true for E,.;. Since | Ey! y, = | Fen + E| Mm 
we get, using (6), 

a(k,) = a(ly1) + a(E) + a( ky 1°), (8) 


while, by the induction hypothesis, 


"_fh- ; 
a(E,y.1) + (-1)""' = 3S (' b ') fa(E*) + (-—1)"~*], 


k=) 
: P ® Ke “. h iy l *kK+1 / = 
a(by1:E) + (-1)""2 = b [a(Bet+!) + (—1)" |. 
k=1 


{Inserting these two identities into (8), we obtain (7), q. e. d. 
It follows from Theorem 5 that there exist a polynomial v(l; IN) in 1 
of degree m and an integer ly (depending on IR) such that 


dim | Fal gp = (i; Mt) — 1, for! = lh. (9) 


In fact: choose a positive integer 4, such that | E,, ~ K| yn is sufficiently 
ample and set 1 =h+. Then, denoting by S an irreducible non- 
singular variety belonging to | /,, — K/| gn, we get, using (5), 


dim | E,| m = dim |K + E,+ S| yp = a(t) + a(F,) + a(S) + a(k,S) — | 


while, as (7) shows, a(/2,) or a(/,S) is a polynomial in h of degree m or 
m — |, respectively. This shows the existence of a polynomial v(/; I) 
satisfying (9). 

This polynomial v(/; YN) is called the postulation of M. The arithmetic 
genus® p,(M) of WM is defined by 


p(M) = (—1)™[v(0O; MY) — 1]. (10) 


LEMMA 2. For sufficiently large |, we have 


v(l; S) = dim | E,| y — dim | E, — S| yy. (11) 


Progf: For large 1, we have v(/; S) = dim |E,S\5 + 1, while, since, 
as the Corollary to Theorem 4 shows, | /;) 4, cuts out on S a complete 
linear system, dim |E,S|s is equal to dim | E,! y, — dim | FE, — S| y, — 1 
Hence we get (11). 

From (11) we get immediately the following 
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THEOREM 6. TSI T\m, then vil; S) = vl; T) and p,(S) = 
p(T). 

THEOREM 7. Assumé that the intersection S-7T is an irreducible non- 
singular variety with the intersection-multiplicity 1. If | U\y = |S +7), 
then we have 

vl; U) = v(l; S) + vl; T) — o(l; S-T), (12) 
pal VU) = palS) + palT) + palS:T) (13) 

Proof: Considering S-7 as a divisor on 7, we get, by (11), o/; ST) = 
dim | £,7) » — dim | /,7’ — ST), for large /, while, by virtue of the 
Corollary to Theorem 4, 

dim | E,7| 7 = dim | E,| gy — dim | E, — T| gq — 1, 
dim | E,T — ST|, = dim | E, — S| m — dim |E, — S— T|y — 1 
for sufficiently large /. Hence we obtain 
— dim E, S|, —dim | FE, — 7} gy, + 
dim E, —_ Ul mM 


v(l; S-T) = dim | F,| yp 


ti 


Using (11), we obtain (12) from this immediately. The formula (13) is 
an immediate consequence of (12). 

3. Let S, be a projective space of dimension d and V4.4, Vax -*-, 
Vax °** be irreducible non-singular subvarieties of S, such that each 
V4. is the section of V4_,,; cut out by a hypersurface &’,, of order n, 
(in particular Vy, = ©’,,). Denoting by ©, a general hypersurface of 
order /, we infer readily that the canonical divisor on Vy_, is given by 
Kaen = €u- Vax, where sy = m + m + +n, —d = 1. Since 


&,-Ve ws nee Va KI Va — = Y, at Ke k+l Vi , _ and v1; Va m1) = 
0 > 6 ? 


dim | @,- Va—x41! V ; + | for large /, we get, by (11), 


d—k + 


o(l; Vez) = v(l, Va-esi) — v(l — ng; Ve—pa). (14) 


d 
v(—l —d — 1; Sa) = (—1)%(; Sy). Starting from this, we get, with 
the help of (14), the formula 


It follows from the well-known formula v(/; Sy) = we that 


o(—l + say Vax) = (—1)4 "0; Vax) (15) 
successively for k = 1, 2,3, .... Now, using (4), we get, for sufficiently 
large /, 

v(l + se; Vax) = dim | G4,- Va_s Pe + ] 
= dim | Ka. + ©: Va ‘ly + | 


d—k 
a a(Va k) + a(&,- Va k)s 
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while, as (7) shows, a(@,- V__,) is a polynomial in / whose constant term 
equals (—1)4*. Hence, using (15), we obtain 


a(Va_-x) + (—1)4* = o( 54; Van) = (—1)4*0(0; Vax). 


This proves that p,(Va_,) = a(Va_,) for the varieties Vg_, = ©'n, C'n,: 
oe ee oe oS eee ® 

4. Now we shall prove the identity p,(V,) = a(V,) by induction on n. 
In the simplest case in which n = |, this identity is an immediate conse- 
quence of the theorem of Riemann-Roch. We assume therefore that the 
identity p,(V,) = a(V,) is true for n < m — | and consider an arbitrary 
but fixed irreducible non-singular variety V’,,, of dimension m 2 2 imbedded 
in a projective space S, of dimension d >2m+ 1. Take d-—m— 1 
general hypersurfaces ©,,, ©, ..., Engom— Of sufficiently high orders 
Ni, No, .» Na-m—1 Containing V», and set Vax = €n,°€,,° 
(k= 1,2,...,d—m-—1). Then, as one readily infers from the prin- 
ciple of counting constants and the theorems of Bertini, each V4_, is an 
irreducible non-singular variety. Furthermore VV _, is the section of 
Va-eyi Cut out by the hypersurface ©,,. Hence, denoting by 6,(Vq_,) 
the v-Betti number of V4_,, we get, by a theorem of S. Lefschetz,’ 


b,( Vin+i) ” b, ( Vin+2) = : = b, ( Va 1) _ b,(Ga), for Vv < m. 


In particular b)(Vingi) = b:(Sa) = 0, be Ving) = b2(Sa) = 1, and conse- 
quently bom( Vingi) = be Vingi) = 1. Set WE = Vinyi and denote by FE, 
the section of Yt cut out by a general hypersurface &, of order h. Again 
set 6S) = p,(S) — a(S) for an arbitrary irreducible non-singular sub- 
variety S of Yt of dimension m. Then, comparing (13) with (6), we get 
6(U) = 6(S) + 8(7), for |Ulm =|S+T7|m, (16) 

provided that | S| y and | U| m are both sufficiently ample and that S-7 
is an irreducible non-singular variety with the intersection-multiplicity 1, 
since p,(S:7) = a(S-7) by the induction hypothesis. Now, consider the 
subvariety V, of P= Vingi as a 2m-cycle on YN. Since be,(M) = 1, we 
can find two positive integers g, r such that qV,, ~ E£,, where ~ denotes 
the homology with respect to integral coefficients, while it follows from 
b,(M) = O that” qV,, ~ EF, implies |qVn|y, = | E,|m. Choose h so large 
that |kV,, + E,|y is sufficiently ample for k = 0, 1, 2, ...,q and let S, 
be a general member of | kV, + E,| a». Then, since | Sel a = | Sy1+ 
Vinl gn, we get, by (16), 6(S,) = 6(S;,1) + 6(V,), and therefore 6(S,) = 
q:6(Vm) + 8(So), while |S.| = | E,yn| yn and |Solg = | Ex|g. Conse- 
quently we obtain qg:6(V,,) = 6(/£,,,) — 6(/,). On the other hand, the 
result of Section 3 shows that 6(/,) = 0 for every E, = ©,,-©,,° 

Eng om 1°G Hence we get 6(V,,) = 0, proving that pa(V,) = a( Vn), 


q.e. d. 
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A GENERALIZATION OF A THEOREM OF HODGE 
By D. C. SPENCER 
PRINCETON UNIVERSITY 
Communicated by S. Lefschetz, April 30, 1952 


Let M be a real Riemannian manifold of dimension n and of class C”. 

At an arbitrary point pe M choose an orthonormal basis ex), k = 1, 2, 
a: k 

.., n, for the tangent vectors at p. Given any real numbers a, k = 1, 


2, ..., m, Satisfying 


the unit vector 


n 
A= Do akey) 
k—-! 


determines a unique geodesic issuing from p. Let r denote distance from 
p along this geodesic, and introduce the normal coordinates x“ = a’r, 
k = 1,2, ..., . These coordinates are obviously determined up to an 
orthogonal transformation. 

Let S,(p) denote the geodesic sphere of radius ¢« and center at the point 
p. We say that M is uniformly smooth if: (i) there is a positive number 
n such that any two points of J whose distance does not exceed n can be 
joined by a unique geodesic; (ii) throughout any S,(p) the partial deriva- 
tives up to the fourth order, with respect to the above coordinates at p, 
of the metric tensor are bounded by a number independent of p. 

We assume that M is uniformly smooth and metrically complete. Then 
the following two theorems are true: 
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(A) Any differential form ¢ of degree g, 0 <q <n, which has finite norm 
and satisfies Ag = 0, also satisfies dg = bg = 0. 

(B) Given any closed q-form a which has a finite norm, there exists a unique 
norm-finile harmonic form which has the same periods on all q-cycles with 
compact carriers. 

Theorem B is a generalization of the classical theorem of Hodge. Let 
C be the group of closed qg-forms ¢ having finite norm (cocycles), and let 
B be the subgroup of C composed of forms g¢ = dy (coboundaries). The 
quotient group // = C/B defines a cohomology group of dimension q 
which is isomorphic to the space of harmonic q-forms. If M is compact 
(closed), // is just the cohomology group of dimension q in the usual sense. 

Theorems A and B are false if the hypothesis of metric completeness is 
dropped. A proof of A and B based on the heat-equation method intro- 
duced by Milgram and Rosenbloom in reference 1 will appear in the 
Annals of Mathematics. 


' Milgram, A. N., and Rosenbloom, P. C., ‘Harmonie Forms and Heat Conduction 
I: Closed Riemannian Manifolds,”’ Proc. Nati. Acapb. Scr., 37, 180-184 (1951). 


POWER-ASSOCIATIVE COMMUTATIVE ALGEBRAS OF DEGREE 
TWO 
By L. A. Koxoris 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF CHICAGOT 
Communicated by A. A. Albert, April 25, 1952 
It has been shown that every central simple power-associative com- 
mutative algebra of degree greater than two and characteristic greater 
than two is a classical Jordan algebra. Such an algebra of degree one 


and characteristic zero is isomorphic to its center.'. In this paper we 


shall show that a simple power-associative commutative algebra of degree 
two and characteristic p > 5 need not be a Jordan algebra. We do this 
by constructing an example. 

Let A = F{x]| be an associative algebra over the field F with unity 
element e and let x” = e and x? = 0. Assume that e is the sum of two 
orthogonal idempotents, u and v. Then A = uA + vA. We take two 
quantities yy and y, and consider the system S = A + uA + wA + yA. 
Further assume that a + ub + yor + yid = 0, for a, b,c, din A, if and only 
iia = b=c=d=0. Since A has dimension p, the vector space S has 
dimension 4p. 

Define multiplication in S by assuming the commutative law, the product 
formulas 
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(yix*) (xu) am [(yx®) xc] = [(yx?)ulx® = '/2(yx*)x* = "/eyx* 1? 
and the formula 


.atpB-l ( 


(yx)(y,x®) = (i — j)(a — B)x : 2) 


where 1, 7 = 0, | and a, 8 = 0,1, ..., p — 1. Relation (1) implies that 
e is the unity element of S. Note that? = 1, a = p,j = B = Oin (2) 
gives (yix?)(yox?) = px’?! Since x? = 0, px’! = 0, and it follows that 
the characteristic of S must divide p. We assume that the characteristic 
ofS is p. 

When Pp is greater than 5, the power-associativity of S is equivalent to 
the identity’ 
4(xy) (sw) + 4(xz)(yw) + 4(xw)(ys) = x[(ys)w + (yw)s + (we)y] + 

y[(xz)w + (xw)e + (wz)x] + 2[(xy)w + (xw)y + (yw)x] + 
w|(xy)z + (xz)y + (yz)x]. (3) 

If the characteristic of S is 3 or 5 other identities must also be satisfied if 
S is to be power-associative. We shall content ourselves with showing 
that S is power-associative for p > 5. This is done by substituting all 
possible combinations of four elements in (3). 

In order to facilitate the proof of the power-associativity of S we make 
some auxiliary computations. For convenience we write ux = g and 


vw = hsothatx = g +h. Then (1) and (2) imply 


[vie )e? (9x9 + xD [(yxDe?] = Velyix®* *)(yx7) + 


V/o(yc*) (yy? 77) = '/2(i — j)(a + B— y)xtTPTr! + 


atBt+y-1 


V/A(a — Dla — B-— y)x atBt+y 1 


= (1 — j)\(a — y)x 
We write the result as 
[(vix™) ge" } (9x7) + (yix™ [Cyr e"] = [ (yee) (yx7) |x’. (4) 


Another relation may be obtained by replacing g with h. Direct computa- 
tion shows that 


(yx) (yy) (ye) + [xe Dye) (yn?) = 


(i — j)(2a — B — y)yxttFt 7! (5) 
and 
[(yix™)g? 127 + [(vex™) 7 Je" + (ye) (927) = yart FOF, (6) 
First consider the case where one of the variables in (3) has the form 
y.x*. If the other variables are powers of x then we use the fact that 
(y.x%)x® = yx? to see that (3) is satisfied. If a power of g or h is used 
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in place of a power of x then (1) implies that each term of (3) is halved so 
we have disposed of this case. 

Let us proceed to substitutions involving exactly two variables of the 
form y,x*. The subscripts of these two elements must differ if we are to 
obtain non-zero products. Because of linearity and of the orthogonality 
of g and h we need only consider the two sets of values x = y.x%, y = yor", 
z= g’,w=h’,andx = yx%,y = yr®, 2 = g’,w = yg’. After making the 
first of these substitutions and using (4) we need to show that 


A{ (yix™)g7}[(yox®)h?] + 4[ (nix )h?][ (yor) 27] = (ix) ({(yox®) gh? + 
[(yox®)h?]27) + (yox®) (L(x 27 Jh® + [Ox Dh lg?) + 
([ (yi) (vox) Jx9) 27 + (L Cyx®) (yor®) |x a’. 
We use the multiplication table of S to reduce this equality to 
(yx *7) (yor? T*) + (yux®t)(yox® T%) = 1/2( yx) (yor? F744) + 
'/o(yox®)(yix®F7*9) + [ (yx) (yor®) [x7 F?. 
Then (2) implies that this expression holds in S. Relation (3) with values 
x = yx", y = yor", s = g’, w = g’ gives 
4[ (yx) (yox”) (272°) + 4[ (nix) 27 1[ (yor) 2"] + 4[Onx™) 2°] (yor®)g7] = 
(yx) (yor? F749) + (yox®)(gix® 749) + [(yx*) (yor*) |g??? 
where we have used (4) and (6). This simplifies to 
da — B)g*tOtrtO-1 4 (qa + y — B— Sxtthtrtent 4 
(a+t6— B- y)xetetrtent =(a-p~y7- Byxetoty tent 
(—BH+ at 7 t d)xttPtr7te-! 4 da — BygrtPtrte-) 


which is satisfied identically. 

If three of the variables of (3) have the form y,x*, we need only make 
the substitution x = y.x%, y = yor", 2 = yor’, w = g’. We apply (4) and 
(5) after making the substitution. Then (3) is equivalent to 


$[(yix™)(yor®) 1[(yox 2°] + 41 (x) (yor) }[(yor®) 2°] = 


(vox”) ([ (yx) (yor?) }(0° + gi) + (vex?) ({ (ix) (yor®) ](x? +g”) + 


g"((2a — B — y)yox" tA+y— J 


Every summand of this expression is a multiple of yox* *°*7*?~! and com- 
putation of the coefficients gives 2(a@ — 8) + 2(a@ — y) = */x(a — y) + 
3/o(a — B) + '/x(2a — B — y) as desired. 

There remains the case with all four variables of (3) having the form 
yx". The only set of values of this type giving non-zero products is 
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x= yx% y = yx, 2 = yor’, w = yor. We then wish to show that 


4[ yx") (yor”) | (yx) (yor?) | + AL (yx) (yor?) | (oe) (yor) ] = 


(yix®)[(2B — y — dyer? F7T9>"] + (yx?) [(2a — y — b)yoxt@ th ™FE>1] + 


(yox") [(a + B — 2)yx™F PFO") + (yor*)[(a + B — Zy)yeettetr—!y, 


= ‘ e . 4 ~§ 
Each term is a multiple of x**°*7*? 


shows that (3) is satisfied. Thus the algebra S is power-associative for 


and a check of the coefficients 


any characteristic p > 5. 

The algebra S is simple. For suppose B is a non-zero ideal of S. If 
u = g’ isin S then 2(y.x)g° = yx and (y.x)(yox"”) = x” = eisin S. There- 
fore Bis S. Similarly, if v = h° isin B, B = S. If x* isin B, it follows that 
[(y.x")x*]yo = ax*~'isin B. After a finite number of applications of this 
argument we have x° = ein B. In the same manner we may show that 
g* in B implies B = S. Whenever y,x* is in B then (y,*)(y") = (@ - 
jax*~'isin B. We take j # j and then have x*~'in B. Then B cannot 
be a proper ideal. 

In the decomposition of S relative to the idempotent u, S,(1) = Fig] = 
uA, S(O) = Flh| = vA, and S,('/2) = yA + yA. It follows from (1) 
that S is u-stable; that is, S,('/2)S,(A) S S,('/2) for’ = 0,1. However, 
Sis not stable. We prove this statement by showing that the subalgebra 
T = A + wA + yA is not stable. The subalgebra 7 has order 3p and 
is obviously simple. Consider the element w = yx + '/2yx° in 7. Then 
w? = (yx)(yor”) = x° = e and [(e + w)/2]? = (e + w)/2. That is, 
z= '/(e + wx + '/syox") is an idempotent. Also, e — z = (e — w)/2 is 
idempotent, e = (e — z) + zand (e — z)z = z — 2? = Oso T has degree 2. 
We have (yix")2 = '/oyix" and (2x)z2 = '/o(x + yix? + '/oyox)e = a(x + 
vix? + '/oyox + yx? + '/oyoxe + x) = ax. Therefore y,x° is in 7,('/2) 
and sx = '/o(x + yx? + '/oyox) is in 7,(1). We compute (y,x°)(2x) = 
'/s(yix — '/ee) and [(yix") (2x) 2 = '/a(yux — '/2e — '/oyie + '/2e — '/eyor”) = 
"44 /oysx — '/gyox®). Since [(yix") (2x) |z # '/2[(y1x") (2x) ], T, and therefore 
S,is not stable. Every Jordan algebra is stable so S and 7 are not Jordan 
algebras.* 

+ This work was carried out in part with the aid of the Office of Naval Research. 

' These results are given in Albert, A. A., Trans. Amer. Math. Soc., 69, 503-527 
(1950) with the restriction that the characteristic be different from 2, 3, and 5. The 
results were extended for characteristic 3 and 5 in the author's unpublished Ph.D. 
dissertation. 

2 See Albert, A. A., Summa Brasiliensis Mathematicae, 2, 21-33 (1948). 

3 The subalgebra 7 was originally constructed independently of S and in a much more 
complicated form. The author is indebted to A. A. Albert for noticing that 7 is a sub 
algebra of S and that S is not stable although u-stable. 
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SOME RESULTS ON ABELIAN GROUPS' 
By IRVING KAPLANSKY 
DEPARTMENT OF MATHEMATICS, THE UNIVERSITY OF CHICAGO 


Communicated by A. A. Albert, April 25, 1952 


1. Introduction.— The author is engaged in the preparation of a mono- 
graph on infinite abelian groups. Besides exposition of known theory, 
it is planned to include some new material. Since publication will pre- 
sumably be somewhat delayed, it seems appropriate to make this brief 


announcement of the new results. 

2. Definitions.-Let G be an abelian group. A subgroup of G is fully 
invariant (characteristic) if it admits every endomorphism (automorphism) 
of G. A subgroup // of G is divisible if nT = H for all n. G is reduced 
if its only divisible subgroup is 0. For most questions it suffices to con- 
sider reduced groups. G is primary (for the prime p) if every element 
has order a power of p. For a reduced primary group G we define a trans- 
finite series as follows: G,,,1 = pG, for every ordinal a, and G, is the 
intersection of the preceding subgroups if a@ is a limit ordinal. If a non- 
zero element x lies in G, but not in G,4;, we write h(x) = a. We write 
h(O) = , with the agreement that © exceeds any ordinal. Let P be 
the subgroup of elements of order p. The dimension of (P A G,)/(P 
G,41), regarded as a vector space over the integers mod p, is the ath 
Ulm invariant of G, written f(a). If G is countable, it is characterized 
by its Ulm invariants. 

3. The Ring of Endomorphisms.—-TueoreM 1. Let G and H be primary 
abelian groups, Ek and F their rings of endomorphisms. Then any 1so- 
morphism of E and F is induced by a group isomorphism of G and H. In 
particular, every automorphism of E ts inner. 

This theorem was proved by Baer,’ for groups of bounded order having 
at least three linearly independent elements of maximum order, with the 
aid of his theory of the lattice of subgroups. The author’s proof uses 
ring-theoretic techniques, reminiscent of the classical treatment of auto- 
morphisms of matrix rings. 

4. Characteristic Subgroups.—In Theorems 2 and 3 we confine ourselves 
to a reduced primary group G, but the theorems are easily extended to the 
non-reduced case. Let a; be a monotone increasing sequence of ordinals; 
it is allowed to be © from some point on. The set of x satisfying h(p'x) 2 
a, for all 7 is manifestly a fully invariant subgroup of G.* 

THrOREM 2. Let G bea reduced primary group such that every two elements 
can be embedded in a countable direct summand (this is in particular true if 
G itself is countable, or if it has no elements of infinite height). Then every 
fully invariant subgroup of G has the form just described. The lattice of 
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fully invariant subgroups is distributive and even satisfies the infinite dis- 
tributive law A MN (UB, = V(A 1 B,); the dual infinite law is, however, 
satisfied if and only if G has no elements of infinite height. 

Characteristic subgroups are fully invariant, with an exceptional case 
generalizing that discovered by Shoda‘ and Baer.° 

THEOREM 3. The hypothesis on G is the same as in Theorem 2. Then 
the following two statements are equivalent: (a) G has a characteristic subgroup 
which is not fully invariant, (b) p = 2, and there exist ordinals a, B with 
B > at] such that the Ulm invariants f(a) and f(B) are both 1. 

5. Complete Modules.—In this final section we are concerned with 
modules over the p-adic integers (any complete discrete valuation ring 
would do in place of the p-adic integers). If / is such a module, 
we say it has no elements of infinite height if the intersection of p'M 
is 0. These submodules may then be taken as neighborhoods of 0 
for a topology, and we say that \/ is complete if it is complete in this 
topology. 

THEOREM 4. A complete module over the p-adic integers is the completion 
of a direct sum of cyclic modules. The cardinal numbers giving the number 
of cyclic summands of each order are a complete set of invariants.® 

From Theorem 4 it is possible to deduce the answer to a question raised 


by the author in reference 7. 
THEOREM 5. Any module M over the p-adic integers has a direct summand 
of rank one. In particular, M is indecomposable if and only tf it has rank 


one. 

Theorem 4 can also be applied to the problem of determining which 
abelian groups can be compact. Let G be a compact abelian group. It is 
known that its component C of the identity is its maximal divisible sub- 
group. Hence C is (algebraically, not topologically!) a direct summand 
of G. It is also known that the totally disconnected group G/C is the 
complete direct product of groups //,, each of which is a complete module 
over the p-adic integers. Thus G belongs to a class of abelian groups 
which is completely characterized by a set of cardinal numbers. The 
author has not fully solved the problem of determining which sets of 
cardinals are eligible to be those of a compact group, but this appears to be 
of lesser interest. Indeed it would seem to be a good idea to call a group 
G “algebraically compact” if it has the form G = C @ //, where C is 
divisible and // is a complete direct sum of complete modules over the 
p-adic integers. Then at any rate negative theorems carry greater force 
when they assert that a group is not algebraically compact. Sample 
theorem: an abelian group with an infinite cyclic direct summand is not 
algebraically compact. 

' This work was supported in part by the Office of Naval Research. 

2 Ann. Math., 44, 192-227 (1943). 
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§ This deseription of fully invariant subgroups is a convenient modification of that 
of Shiffman, Duke Math. J., 6, 579-597 (1940) 

* Math. Zeit., 31, 611-624 (1930) 

® Proc. London Math. Soc., 39, 481-514 (1935) 

* The torsion-free case of Theorem 4 was discovered by Isidore B. Fleischer and 
appears in his dissertation. 

7 Trans. Am. Math. Soc., 72, 327-340 (1952). 


NATURE OF INHERITED RESISTANCE TO VIRUSES 
AFFECTING THE NERVOUS SYSTEM 
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THe CHILDREN’S HospiraL RESEARCH FOUNDATION, UNIVERSITY OF CINCINNATI 
COLLEGE OF MEDICINE, CINCINNATI, OHIO 


Read before the Academy, April 28, 1952 


Introduction.—The vast majority of viruses which attack the nervous 
system of human beings and animals produce recognizable disease and 
death in only a small proportion of infected individuals. A variety of 
factors pertaining to the virus, the host, or both, may influence the out- 
come of infection. The studies to be presented here were designed to 
elucidate the intricate biological mechanisms which form the basis of this 
important phenomenon. 

Previous studies by other investigators have established that plants and 
animals may possess an inherited resistance to various infectious and 
noxious agents, including viruses. The studies of Lynch and Hughes! 
with the virus of yellow fever and those of Webster? with the viruses of 


louping ill and St. Louis encephalitis provided the first experimental evi- 


dence that the genetic constitution of the host can determine the outcome 
of mammalian viral infections. A most important ‘contribution to this 
subject was Webster's demonstration that the inherent resistance or sus- 
ceptibility of mice to the virus of St. Louis encephalitis was correlated with 
the level of viral multiplication in the brain, not only in the intact animal? 
but also in simple cultures containing the minced brain tissue.‘ Earlier 
attempts to establish the manner in which such resistance is inherited by 
animals, and to segregate the factors which are genetically affected, were 
complicated by the fact that no uniformly resistant animals were available. 

Experimental Results.—In 1944, the author accidentally discovered that 
the albino mice which had been bred for over 25 years at the Rockefeller In- 
stitute at Princeton, N. J., were 100 per cent resistant to the 17 D strain of 
yellow fever virus that is widely used for human vaccination. Swiss mice, 
intracerebrally inoculated with this virus, invariably die after exhibiting 
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paralysis and other signs of involvement of the nervous system. In the 
Princeton Rockefeller Institute mice (henceforth called PRI), intracerebral 
inoculation of the 17 D virus, in the largest as well as the smallest doses, 
results in multiplication at a level only '/ 0.000 to '/ 00,00 of that achieved 
in Swiss mice, and in the development of specific antibodies—but all mice 
survive without showing any clinical signs of disease. The possibility 
was considered that a latent virus in the PRI mice might interfere with 
the multiplication of the yellow fever virus, but none was found in tests 
with brain, viscera, intestines and feces. This circumstance presented 
an unusual opportunity for genetic analysis by the classical methods of 





MECHANISM OF INHERITANCE OF RESISTANCE TO 17D YELLOW FEVER VIRUS IN MICE 





(100% die) SUSCEPTIBLE RESISTANT (100 % remain well) 
SWISS (300) PRI(OO) 





Fs 


BACKCROSS (79) 





Letters within circles=genetic 
formula fitting observed data 


(100) No. of mice tested 

















FIGURE | 


interbreeding. Such studies were begun in 1950, and this communication 
is a progress report of the results obtained thus far. 

Figure 1 presents a diagrammatic summary of the tests with the 17 D 
virus. The sex of the parents and progeny is not indicated because no 
evidence of sex-linkage was found. Since all F; progeny survived without 
signs of disease, it is clear that resistance to this virus is inherited as a 
dominant. The approximate 3:1 ratio of resistants to susceptibles in the 


F, progeny, the exact 1:1 ratio in the backcross of F; to susceptible Swiss 
mice, and the 100 per cent resistance in the backcross of F, to resistant 
PRI mice, are all in accord with the Mendelian laws for a single pair of 
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autosomal genes. However, in 1934, Sewall Wright® demonstrated that 
such ratios are not critical criteria for the operation of a single pair of genes 
and that the critical experiment must test the genetic nature of the back- 
cross and F, progeny by breeding them with the recessive stock. This was 
done, as shown in figure 1, and the results obtained were again compatible 
with a single factor hypothesis. It is realized, however, that it is difficult 
to rule out other possibilities involving multiple genes, which do not assort 
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FIGURE 2 


This figure shows the mortality resulting from tests on parental and F, groups of mice. 
The genetic aspects are discussed in the text. 


at random in the gametes but are linked with only a limited amount of 
crossing over. Thus, if resistance to 17 D were dependent on two pairs of 
dominant, linked genes, with LO per cent crossing over, the backcross and F, 
progeny and the progeny resulting from crossing these with the recessive 
stock would be modified so slightly from the single factor expectancy that 
only exceptionally large numbers of animals might yield results of sufficient 
statistical significance to permit a decision. However, whatever the exact 
number of genes may be, it is clear that the resistance of PRI mice to the 
17 D virus depends on a factor, which depresses viral multiplication, and 
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is inherited as an “‘all or none’”’ character, because the progeny which segre- 
gated out of the various crossings were phenotypically like the parent 
resistant or susceptible stock, i.e., the resistant mice showed no signs of 
disease and the susceptible ones exhibited the same incubation period and 
high level of viral multiplication as the Swiss mice. 

All this held good as long as the tests were restricted to mature animals. 
When newborn, |-day or 2-day old PRI mice were inoculated with 17 D 
virus, all mice died after an incubation period of 10 days or longer. PRI mice 
inoculated at 3 to 5 days of age behaved irregularly, and beyond 5 days of 
age, the suckling mice behaved like the mature animals. Tests for viral 
multiplication in the succumbing suckling PRI mice showed that the peak 
levels were as low as in adult PRI mice which remained well. It was ap- 
parent, therefore, that at low levels of viral multiplication the tissues of the 
newborn mice were more vulnerable than those of older mice. 

The concept of cellular vulnerability again came under consideration 
when it was found that the ‘French neurotropic’’ strain of yellow fever 
virus killed a certain proportion of adult PRI mice. This strain of virus 
has had over 250 serial passages in the brains of Swiss mice, and although 
its level of multiplication is not significantly higher than that of the 17 D 
strain, it kills mice more rapidly. The results summarized in figure 2 
show that 24 per cent of adult PRI mice died after intracerebral inoculation 
of the ‘French neurotropic” virus. The virus multiplied at the same low 
level in the PRI mice that died as in those that remained well, and 
both of these levels were about the same as the peak levels of the 17 D 
strain. There was no evidence that the virus particles which grew out in 
the brains of the dying mice were different from those in the brains of mice 
which remained well, since passage to new PRI mice yielded a mortality 
of only 20 per cent. It appeared, therefore, that some of the PRI mice 
were more vulnerable than others at similar low levels of viral multipli- 
cation, and the question arose as to whether cellular vulnerability might 
also be genetically determined. 

Since the multiplication-depressing factor is inherited as a dominant, 
it is possible to compare the cellular vulnerability of the F; progeny (Swiss 
X PRI) mice with that of PRI mice, without having it complicated by 
the factor of viral multiplication. In repeated tests, in which approxi- 
mately one million infective doses of “French neurotropic’’ virus were 
inoculated intracerebrally in groups of 20 to 30 mature PRI mice, the 
mortality ranged from 20 to 30 per cent, with an average of 24 per cent. 
Crosses made between PRI and the highly inbred CFW, Swiss mice pro- 
duced 84 F, animals, which exhibited an average mortality of 70 per cent 
after intracerebral inoculation of the same large dose of “French neuro- 


tropic” virus. In litters of 7 or more mice the mortality was never less 
than 50 per cent, and in 2 families of 8 mice, each resulting from a cross 
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of PRI females with Swiss males, the mortality was 100 per cent. In 
general, it made no difference whether the mother was of PRI or Swiss 
stock. This result, which is highly significant statistically, indicates that 
the character of the parents influences the behavior of the progeny. 
When F, mice were backcrossed to PRI mice (not shown in figure 2), 
60 per cent of the progeny (12/20) died following inoculation of the 
French neurotropic virus. These results strongly suggest that cellular 
vulnerability to the ‘French neurotropic’ strain of yellow fever is also 
genetically determined. The fact that PRI parents, which were proved 
to be resistant, i.e., possessed a low cellular vulnerability, yielded 18 per 
cent of highly vulnerable progeny suggests not only that high cellular 


TABLE 1 
SELECTIVE ACTION OF MULTIPLICATION INHIBITION FACTOR OF PRI MICE ON VARIOUS 
VIRUSES PROLIFERATING IN MOUSE BRAIN 
ACTION VIRUS 
| Yellow Fever 
Dengue Fever 
Multiplication inhibited. PRI | West Nile Fever 

mice completely or partly re Japanese B Encephalitis 

sistant St. Louis Encephalitis 
Russian Spring-Summer Encephalitis 
| Western Equine Encephalitis 
Eastern Equine Encephalitis 
Venezuelan Equine Encephalitis 
Poliomyelitis 

Multiplication not affected. | Mouse Encephalomyelitis —‘‘TO” 

mice fully susceptible Rabies 
| Lymphocytic Choriomeningitis 
| Herpes Simplex 
| Vesicular Stomatitis 
Rift Valley Fever 


vulnerability might involve a recessive gene or genes but also that many 
of the unselected PRI stock are heterozygous for this character. The F, 
progeny resulting from crossing resistant PRI mice with Swiss mice yielded 
46 per cent of highly vulnerable animals, which again would suggest that 
the Swiss mice are probably homozygous for the recessive genes, responsible 
for high vulnerability, and that the resistant (low vulnerable) PRI mice 
are predominantly heterozygous. The progeny resulting from crossing 
susceptible PRI with resistant PRI contained a lower percentage (31 per 
cent) of highly vulnerable individuals than in the tests in which Swiss 
mice were used; the difference between this result and the 18 per cent of 
vulnerable progeny produced by resistant PRI parents is in the right 
direction but not large enough for statistical significance. A great deal 
more work will have to be done on segregating the unselected PRI mice 
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by individual matings with Swiss mice before further genetic analysis of 
the factor of cellular vulnerability will become possible. While the viral- 
multiplication factor yielded ‘all or none’’ results, there is an indication 
that cellular vulnerability is not an ‘‘all or none” character, since a small 
proportion of mice develop weakness or extensive paralysis of the extremities 
but do not die. 

The genetic factor in PRI mice, which depresses multiplication of the 
yellow fever virus, was found to exert a similar effect on the viruses of 
dengue fever, West Nile fever, Japanese B, St. Louis and Russian spring- 
summer encephalitis, but was without effect on a large group of other 
viruses listed in table |. It 1s of interest to note that the viruses, which 
are affected by this genetic factor, are also linked together by a chain of 
common antigens,® not shared by the other viruses. 


TABLE 2 
OCCURRENCE OF VARIANT STRAINS OF VIRUS CAPABLE OF OVERCOMING INHERITED 
RESISTANCE OF Host 
PERAK 
BEHAVIOR IN PRI MICE MULTIPLICA 
NO. OF PERAK TION IN 


PASSAGES IN MORTALITY MULTIPLICATION BRAINS OF 
VIRUS STRAIN SWISS MICE o// IN BRAIN SWISS MICE 


Japanese B-— Korea 13 2.4.3.0.3:9 8 0" 
encephalitis Nakayama : 50 a 8 


Nakayama 100 7.6 95 
St. Louis Winkler s 6 & 
encephalitis Webster No.3) ‘‘Hundreds’’ over 90) , 9 0) 


17 years 


* Reciprocal of log of LDyo 


Only the dengue viruses behaved like the 17 D yellow fever virus in that 
all the PRI mice were uniformly resistant regardless of dosage. The 
behavior of the West Nile, Japanese B and St. Louis viruses varied mark- 
edly with the strain of virus and the number of passages it had had in 
Swiss mice prior to test in PRI mice. Some of the strains killed no PRI 
mice, others a varying proportion, and still others killed all the animals. 
Table 2 presents some of the data which indicate that viruses may over 
come the inherited resistance of the host in at least 2 different ways: (1) by 


possessing or developing variants which can kill at low levels of viral 


multiplication, and (2) by the development of variants which are no longer 
inhibited by the inherited multiplication-depressing factor. We cannot be 
certain that these variations in the populations of particles belonging to 
any one virus may not occur in nature, and, the genetically resistant mouse 
can be used to investigate this possibility. For example, it was found that 
2 of the 3 strains of West Nile virus recovered in Dr. Paul's laboratory 
from Egyptian sera,’ killed only a s nall proportion of PRI mice, while the 
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third strain (Egypt 101), in the second mouse passage, killed all PRI mice 
although it multiplied at a level that was only '/ 4099 of that achieved in the 
brains of Swiss mice. 

Summary and Discussion.These studies, in brief, have shown that 
an analysis of inherited resistance to viruses must take into consideration 
not only factors which control viral multiplication but also the phenomenon 
of cellular vulnerability. The available data suggest that cellular vulner- 
ability also may be a genetically controlled character. With some strains 
of virus only one of these factors comes into play, with others, both. 
The demonstration that any one virus does not consist of a homogeneous 
population of particles, and that certain strains can possess or develop 
the capacity to overcome one or the other, or both, of these host barriers, 
poses a problem of considerable interest with regard to the potential peri 
odie emergence of viruses with special epidemic or epizootic properties. 
These data provide a model which indicates how the viruses responsible 
for encephalitis and poliomyelitis may behave in populations of mixed 
or highly selected genetic constitution. The data also provide an explana- 
tion for the occurrence in some of the highly endemic areas of Africa of 
population groups which are almost uniformly resistant to yellow fever. 


The mixed occurrence of severe and mild infections among the South 
American Indians is in accord with the hypothesis that yellow fever was 
imported there from Africa after the discovery of America, and that not 


enough time has elapsed for the virus to kill off a significant number of 
those carrying the recessive genes for susceptibility. Finally, these studies 
provide us with a model of how nature prevents viruses from causing disease 
by allowing them to multiply just enough to produce immunity, but not 
enough to damage the tissues—a model that we may well investigate in 
great detail because an understanding of the underlying mechanisms would 
provide us with the ideal procedures for the chemoprophylaxis of virus 
diseases. 
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“The underlying physical laws necessary for the 

mathematical theory of a large part of physics and 

the whole of chemistry are thus completely known, 

and the difficulty is only that the exact application 

of these laws leads to equations much too compli- 
cated to be soluble.’’—Dirac* 


In April of 1950, at the Detroit ACS meeting, a group of chemists and 
physicists fell into a discussion of the inadequacies of valence theory, and 
what might be done to improve it. The need felt most acutely was a 
reliable table of the diflicuit integrals which inevitably turn up in valence 
calculations. Everyone agreed that a conference on such integrals would 
be desirable, both to find out which integrals would be most useful, and to 
avoid duplication of effort in computing them. Actually, the evaluation 
of integrals is not the only difficulty in quantum-mechanical valance theory, 
and as the idea of a conference developed, the range of topics was broadened. 
R. S. Mulliken of the University of Chicago took this idea and carried out 
the considerable work of organizing a conference. In this he was aided by 
the advice and cooperation of D. A. MacInnes, Chairman of the National 
Academy of Sciences Committee on Scientific Conferences. And so, with 
the sponsorship of the National Academy of Sciences and the financial 
support and encouragement of the Office of Naval Research, a Conference 
on Quantum-Mechanical Methods in Valence Theory was held on Septem- 
ber 8 to 10, 1951, at Ram’s Head Inn, Shelter Island, Long Island, with 
Mulliken as Conference Chairman. 

Invitations were issued to a number of chemists and physicists, both in 
the United States and overseas. The following attended: 

M. P. Barnett (Rayon Research Associa P. O. Léwdin (Uppsala) 

tion, London) ID. A. MaclInnes (Rockefeller Institute, 

T. H. Berlin (Johns Hopkins) New York) 

B. L. Crawford, Jr. (Minnesota ) H. Margenau (Yale) 

C. A. Coulson (King’s College, London ) J. E. Mayer (Chicago) 

H. Eyring (Utah) W. Moffitt (British Rubber Producers’ 
J. O. Hirschfelder (Wisconsin ) Research Association, Herts) 

G. E. Kimball (Columbia ) R. S. Mulliken (Chicago) 

M. Kotani (Tokyo) R. G. Parr (Carnegie Institute of Tech 
J. E. Lennard-Jones (Cambridge) nology ) 
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kK. $. Pitzer (California ) Technology) 

C.C. J. Roothaan (Chicago) L. E. Sutton (Oxford ) 

kK. Riidenberg (Chicago) C. W. Ufford (Pennsylvania ) 
H. Shull (lowa State) J. H. Van Vleck (Harvard ) 
J. C. Slater (Massachusetts Institute of © G. W. Wheland (Chicago) 


The salubrious influence of the setting and accommodations engendered 
a corresponding freedom and informality in the discussions. These in- 
cluded reports on work completed, summaries of work in progress, indica- 
tions of work contemplated, explanations of general viewpoints and argu- 
ments concerning them, helpful hints to quantum calculators, and concerted 
discussion and planning, particularly in the matter of integral evaluation. 
The diverse items, as collected in the ONR Report on the Conference, ! 
give a picture of valence theory as it stands today. 

The contributions to the discussions and the items which have been col- 
lected in the ONR Report seem to fall into six groups. The first group 
contained essentially straightforward applications of present-day theory. 
Most of these dealt with molecular problems; the atom was, however, rep- 
resented by Ufford’s calculations on d° and d°s configurations of Mn IT and 
Fe III, with configuration interaction, The molecules discussed ranged 
from hydrogen to pyrene, in order of increasipg complexity of molecule and 
increasing empiricism of approach. Kotani is studying the polarizability 
of hydrogen, using an accurate wave function of the James-Coolidge type; 
his group is also treating acetylene as a pi-electron system by the LCAO 
molecular orbital method including configuration interaction. The nitro- 
gen molecule is being studied in detail by several groups, including those 
of Herzfeld and Kopineck. 

Shull presented a paper on the charge distribution in the carbon-oxygen 
bond of carbon dioxide which triggered an active discussion.” Several 
workers are now engaged in accurate calculations of dipole moments; use- 
ful comparisons can be made between theoretical and observed values of 


dipole moments and their dependence on internuclear distance, but the 


calculation must be detailed and the comparison critical. 

There was little discussion at the Conference on aromatic molecules, 
which probably reflects a certain satisfaction with present pi-electron 
theory. Shull reported that varying the shielding constant in the 
Goeppert-Mayer-Sklar treatment of benzene does not alter the order of the 
computed electronic levels. Coulson described some calculations on the 
twisting of diphenyl which led to the conclusion that the equilibrium posi- 
tion might be as much as twenty degrees away from the planar. Wheland 
discussed magnetic anisotropy in a number of aromatic molecules, as 
computed by an empirical LCAO method, including overlap.’ 

The second group of contributions considered the quantum interpretation 
of chemical valence concepts—electron pairs, bond energies, bond orders, 
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hybridization, and chemical reactivity. Lennard-Jones indicated how tre 
motion of electrons in molecules can be interpreted in terms of electron pairs 
whose positions are correlated, and showed that a useful approach to molec- 
Berlin pointed 


’ 


ular problems can be based on ‘“‘correlation functions.’ 
out that the space around the nuclei in a diatomic molecule may be pre- 
cisely divided into binding and antibinding regions, presence of an electron 
in the binding region tending to pull the nuclei together, in the antibinding 
region tending to make the nuclei fly apart. Mulliken disclosed the 
latest improvements in his “‘magie formula” for the semempirical calcula- 
tion of bond energies. He emphasized that pi-bonding contributions to 
the bond energies of multiple bonds are probably much larger than gener- 
ally believed, and that isovalent hybridization (hybridization without 
change in valence) and non-bonded inner shell exchange repulsions are often 
important factors in determining molecular stability. Coulson and 
Mulliken discussed the concept of bond order, particularly with respect to 
the difficulties which arise when one attempts to define bond order for 
heteronuclear bonds. ‘The concepts of bond order and hybridization both 
run into complications due to configuration interaction, as was pointed out 
by Daudel for the first and Mulliken for the second. 

Coulson displayed maps of electron density in benzene, based in part on a 
Fermi-Thomas treatment of this molecule; these showed that the pi- 
electrons, on which so much attention has been focused, are actually ‘‘to- 


tally immersed in a sea of sigma-electrons.”* This overwhelming density 


of sigma-electrons must produce overwhelming effects, and the question 
remains why the pi-electron theory of aromatic behavior is so successful. 

Eyring discussed the possibility of calculating by quantum mechanics 
the structures of activated complexes, pointing out that this offered one 
of the greatest rewards of a successful valence theory. With the greater 
reward, the problem carries far greater difficulties, arising chiefly from the 
dearth of a priort knowledge of the equilibrium configuration of the acti- 
vated complex. Even the simplest complex, H,, has not yet been ade 
quately treated. 

The third group of contributions dealt with forces between molecules and 
between non-bonded atoms. First there was a review by Hirschfelder, 
which came to the unhappy conclusion that van der Waals’ forces are still 
better determined from empirical data on gases than from theory. 
He stressed that a new attack on intermolecular forces, using new values 
of the appropriate integrals, would be in order, beginning with the simplest 
case of the inert gases. Next Margenau told of some recent work, including 
a treatment of the London forces between helium atoms which leads to a 
convergent rather than a divergent (asymptotic) series for the potential 
energy as a function of interatomic distance and which gives excellent 
agreement with experiment. 
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Very long range forces, operating over distances of 50 A. or more, have 


recently been postulated by some biological workers. Coulson described 
an investigation of a possible origin of such forces, namely dispersion 
forces between extended oscillators. Deviations from the usual inverse 
sixth-order law were found, but not sufficiently startling ones. The model 
was crude, and the work 1s being extended; the conclusions are tentative. 
A specific type of intermolecular force, distinct from dispersion forces, and 
which he calls charge-transfer forces, was discussed by Mulliken, who 
attributes to these the formation of loose molecular complexes, and the 
interaction of electron donors and acceptors in general. 

Van Vleck discussed a new source of information on intermolecular 
forces, line-breadth in microwave spectra. Line breadths furnish a meas- 
ure of dipole-dipole interactions in the self-broadening of ammonia and 
hydrogen eyanide; of dipole-quadrupole interactions in the broadening of 
ammonia by carbon dioxide; and of dipole-quadrupole and even quadru- 
pole-quadrupole interactions in the broadening of ammonia by nitrogen 
and oxygen, 

The fourth group of contributions was concerned with the problem usually 
spoken of as that of configuration interaction or correlation energy. In 
attacking molecular problems the valence-bond or molecular orbital ap- 
proaches have usually been used. Serious error is inherent in the first ap- 
proximation of either approach; the error in the energy is known as “‘cor- 
relation energy” in the case of the molecular orbital approximation for the 
normal state of most molecules. In either approach, the wave function is 
built up from orbitals (one-electron wave functions), which may be localized 
or non-localized, atomic or nonatomic. No one has yet exhibited a quan- 
titatively satisfactory first approximation. In principle, having chosen 
any approach and enough orbitals of any appropriate type, one can always 
converge to the correct answer by successive approximations. In practice, 
the convergence must be rapid for the method to be useful. For a given 
problem, this depends critically on the approach chosen and the starting 
orbitals used. 

This choice of starting approach and starting wave functions is a core 
problem of quantum-mechanical valence theory. Much attention at the 
Conference was focused on this question. Slater discussed it not only in 
connection with magnetic studies, but also in the general case. Some of 
those present discussed the relative importance of correction terms added 
to conventional first approximations; these included Parr on helium,® 
Mulliken® and Lennard-Jones’ on hydrogen, and Kotani on water and the 
methyl radical. Others discussed the choice of best first approximations 
of the orbital type. Kotani and Eyring each told of investigating new 
orbitals for hydrogen, and Eyring proposed the use of a new sort of ‘‘semi- 
localized orbital’ which he found to work well for both the hydrogen and 
hydrogen fluoride molecules.* 
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Mofhtt, after an incisive critique of the usual approaches, outlined a 
method to avoid the pitfalls of the orbital approximations. Rather than 
computing total molecular energies relative to separated nuclei and elec- 
trons, he computes directly the chemical dissociation energies, that is, 
molecular energies relative to the separated atoms. Instead of starting 
from orbitals, he starts from total atomic wave functions; these are of course 
not known, but by an elegant method, he obtains the molecular energies 
in terms of atomic term values and relatively small perturbation terms. 
Since the perturbation terms are small, they may be calculated from ap 
proximate wave functions, such as orbital functions. The needed atomic 
term values may be obtained either from separate atemic calculations, or 
(a more likely source) from spectroscopic observations on atoms. In 
principle, use of enough atomic term values ensures convergence; these 
must include ionic as well as neutral term values. The question which can 
only be settled by calculation is whether so many term values will be re- 
quired that the method will be unwieldy. Some reported first results on 
oxygen and ethylene are encouraging. 

A fifth group of contributions dealt with mathematical developments. 
Maver outlined a formulation of the Hartree-Fock equations employing 
the charge density matrix and Lowdin reviewed numerical integration 
schemes for these equations. Slater and Ufford reported that IBM 
machines are now in use for the calculation of atomic self-consistent fields, 
but Slater cautioned that direct attack on problems of electronic structure 
is not yet within the reach of automatic computing machines; the machines 
can do much complex arithmetic, but the problem of the molecule has not 
yet been reduced to an arithmetic level. 

The sixth group of contributions dealt explicitly with integrals. In 
clusion of formerly neglected ‘‘small”’ integrals in molecular calculations of 
recent years has led to the realization that “completely theoretical” (as 
contrasted to “empirical’’) valence-theory calculations are capable of 
giving useful results, but with this realization comes an increased need for 
comprehensive table of integrals—in any valid molecular calculation all 
integrals must be included, unto the fifth center and the farthest neighbor. 


If one uses an approach based on atomic orbitals, in any problem one is 


ultimately led to integrals involving these orbitals alone. The integrals 
are of several types, which may be thought of as generalizations of the inte 
grals met in the classical Heitler-London theory of the hydrogen molecule. 
Anyone who has worked through that theory, including the evaluation of 
the exchange integral by Sugiura, will appreciate the task now confronting 
the quantum mechanic. Small wonder that much discussion at the Con- 
ference centered on the topic of integrals, and that much interest was 
aroused by the possibility of cooperative work on the job ahead! 

Groups now working on integrals include those at Cambridge University, 
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where Boys is systematically setting up to compute two-, three-, and four- 


center integrals using automatic computing machinery; the Max-Planck 
Institute in Gottingen, where Kopineck has completed a detailed compila- 
tion of known values for homopolar two-quantum two-center integrals and 
is extending the available values with new computations; Catholic Univer- 
sity, where the program for evaluation of two-center integrals initiated by 
Sklar is being continued with the aid of IBM machines; University of 
Chicago, where Roothaan and Riidenberg, with the cooperation of Shull 
and IBM of lowa State, have well underway a program for computing and 
tabulating two-center integrals by explicit integration; King’s College 
London, where Coulson and Barnett have laid out and begun computations 
with a comprehensive program making use of a scheme they have developed 
whereby a given integral (two-, three-, or four-center) may be calculated 
from a basic set of unit molecular integrals which they propose to tabulate ; 
University of Tokyo, where Kotani is systematically extending his already 
well-known tables of two-center integrals; and University of Uppsala, 
where Lowdin and Lundquist have demonstrated the feasibility of evalua- 
tion by direct numerical integration of each and every integral encountered 
in a given problem. 

More exploratory work needs to be done before a best method for coin- 
puting molecular integrals, if one exists, will be apparent. Probably one 
method will turn out to be best for one type of integral, another best for 
another type. The Conference allowed the several groups of workers 
computing integrals to become acquainted with each other’s work. As a 
consequence, efforts will likely no longer be expended on fruitless forcing 
of a method in a situation to which it is not well suited. The Conference 
also took a step which should greatly help those needing integrals in the 
immediate future: an informal Integrals Committee was established cen- 
tered at the University of Chicago to collect and dispense integral infor- 
mation. 

There is one question that Boys, Kimball and others have recently 
raised: Are Slater atomic orbitals really the best orbitals to use for molec- 
ular calculations? More attention might well be given this question before 
large-scale integral calculations are carried very far; the answer will have 
important bearing on the whole future course of quantum-mechanical 
methods in valence theory. 

Dirac’s statement quoted at the beginning of this summary has long 
been the hope and the despair of valence theoreticians: hopeful because the 
problem is solved “in principle’; despairing because the equations are 
“much too complicated to be soluble.” The encouraging note of the Con- 
ference is that a frontal attack is at last being made on these equations—a 
hopeful, rather than a desperate, attack. 
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* Dirac, P. A. M., Proc. Roy. Soc., A123, 714 (1929). 

' The ONR Report on the Conference consists of numerous informal summaries and 
discussion remarks presented at the Conference, plus some longer papers, a total of 
some fifty-three items. A few of these have appeared in these PROCEEDINGS, 38, 149 
166, 364-378, 400-408 (1952). Copies of the ONR Report may be obtained, as long 
as the supply lasts, from L. M. McKenzie, Head, Physics Branch, Office of Naval 
Research, Department of the Navy, Washington 25, D. C 

* Shull, H., these PROCEEDINGS, 38, 400-408 (1952). 

* Wheland, G. W., and Matlow, H. L., /bid., 38, 364-371 (1952) 

* Coulson, C, A., March, N. H., and Altmann, S., /bid., 38, 372-378 (1952). 

® Taylor, G. R., and Parr, R. G., /bid., 38, 154-160 (1952). 

® Mulliken, R. S., /bid., 38, 160-166 (1952). 

7 Lennard-Jones, J. E., /bid., 38, 496-502 (1952). 

* Mueller, C. P., and Eyring, H., /bid., 38, 149-154 (1952) 
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